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CN ■ Abstract 
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ly-^ I In this paper, we develop a weak-coupling treatment of nonperturbative QCD 

^ . 

I to heavy hadrons on the light-front. First, we present a derivation of quark 

I , confining interaction in light-front QCD for heavy quark systems, based on 
Qh. 
(D . 

r~| , the recently developed light-front similarity renormalization group approach 

. ^ , and the light-front heavy quark effective theory. The resulting effective light- 

X, 

^ , front QCD Hamiltonian Hx at a low-energy cutoff A manifests the coexistence 

of a confining potential and a Coulomb potential. A clear light-front picture 
of quark confinement emerges. Using this low energy QCD Hamiltonian Hx, 
we study heavy hadron bound state equations in the framework of a recently 
proposed possible weak-coupling treatment of non-perturbative QCD. Light- 
front heavy hadron bound states with definite spin and parity are constructed 
and the general structure of the corresponding wavefunctions is explored. A 
Gaussian-type wavefunction ansatz is used to solve the light-front quarkonium 
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bound state equation. We find that the effective coupling constant determined 
from the quarkonium bound state equation can be arbitrarily small so that 
the weak-coupling treatment to heavy hadron bound states in light-front QCD 
is explicitly achieved. Finally, the scale dependence of the effective coupling 
constant is analytically calculated and the similarity renormalization group 
(3 function is determined, from which the running coupling constant in small 
momentum transfer is given qualitatively by a(Q^) ~ . Such a running 
coupling constant is the basic assumption in the successful Richardson QQ 
potential that ensures the existence of a linear confining potential at large 
distance, but now can be obtained from light-front QCD. 
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I. INTRODUCTION 



There are two fundamental problems in QCD for hadronic physics, the quark confine- 
ment and the spontaneous breaking of chiral symmetry. These two problems are the basis 
for solving the low-energy hadronic bound states from QCD but none of them has been 
completely understood. Recently, Wilson and his collaborators proposed an approach to de- 
termine the low-energy bound states in nonperturbative QCD as a weak-coupling problem 
||T|]. The key to eliminating necessarily nonperturbative effects is to construct a low-energy 
QCD Hamiltonian in which quarks and gluons have nonzero constituent masses rather than 
the zero masses of the current picture. The use of constituent masses cuts off the growth 
of the running coupling constant and makes it conceivable that the running coupling never 
leaves the perturbative domain. The weak-coupling approach potentially reconciles the 
simplicity of the constituent quark model with the complexities of QCD. The penalty for 
achieving this weak-coupling picture is the necessity of formulating the problem in light- 
front coordinates and of dealing with the complexities of renormalization. To handle the 
complexities of light-front renormalization, a new renormalization approach, so-called the 
similarity renormalization group scheme, has also recently been developed IHJ^. 

Based on the idea of the light-front similarity renormalization group scheme and the con- 
cept of coupling coherence ||^, Perry has shown that upon a calculation to the second order, 
there exists a logarithmic confining potential in the resulting light-front QCD effective Hamil- 
tonian 0]. This is a crucial finding for a practical realization of the weak-coupling treat- 
ment to nonperturbative QCD. However, the general strategy of solving hadrons through 
the weak-coupling treatment scheme is far to be completed. In this paper, we shall use 
the similarity renormalization group approach to derive a low energy heavy quark QCD 
Hamiltonian, based on the light-front heavy quark effective theory we developed recently 
0J^. The resulting Hamiltonian exhibits explicitly the confining and Coulomb interactions. 
We thereby adopt the idea of the weak-coupling approach to study the strongly interacting 
heavy hadronic structure on the light-front. From this investigation we explicitly provide a 
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realization of the weak-coupling treatment to nonperturbative QCD. 

The reason of choosing heavy hadron systems as a starting example is clear: For light 
quark systems, both quark confinement and spontaneous chiral symmetry breaking play 
an essential role to the quark dynamics in hadrons. However, the success of the chiral 
symmetry description of the low energy hadron physics naturally indicates a chiral symme- 
try breaking scale {A^sb ~ 1 GeV) which is relatively larger than the confinement scale 
(Aqcd ~ 200 MeV). In other words, for light quark systems, it seems to be necessary to 
understand the underlying mechanism of chiral symmetry breaking before we can further 
explore the mechanism of quark confinement. Of course, for the best description, both 
problems should be solved simultaneously in the same picture, but at the moment, this 
will certainly complicate the study on confinement. It would be nice if we could separately 
deal with these two most difficult but fundamental problems in QCD. For heavy quark sys- 
tems, chiral symmetry is explicitly broken so that confinement is the sole nontrivial feature 
influencing heavy quark dynamics. 

Conflnement interaction may not be very sensitive but it is important in describing the 
QCD dynamics of quarkonium spectroscopy and their decay processes. And it should play a 
more important role in heavy-light quark systems, such as B and D mesons. One may argue 
that the mass scales for heavy and light quark systems are different. The heavy quark energy 
cannot run down to the usual hadronic scale of light quark systems due to heavy quark mass. 
Meanwhile, conflning interactions must be energy scale dependent. Apparently, conflnement 
for heavy quark systems could be very different from light quark systems. However, despite 
heavy and light quarks, conflnement arises only from low energy gluon interactions. In 
other words, the conflnement mechanism should be the same for both heavy and light quark 
systems. We thus choose heavy hadron systems without any loss of generality. In order 
to avoid the possible confusion about the different mass scales and to correctly extract the 
conflning interactions in low energy heavy quark dynamics, it is convenient to work with 
heavy quark effective theory (HQET). The HQET is a theory of QCD in l/mq expansion 
0, where tuq is the heavy quark mass. In HQET, the low-energy dynamics is determined 
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through the interacting gluons and heavy quarks by exchanging a small residual momentum 
of heavy quarks, which is of order Aqcd- As a result, within HQET we can indeed explore the 
low energy QCD dynamics for heavy quark systems in the same scale as that for light quark 
systems. Meanwhile, the extension of the study to light quark systems becomes obviously 
straightforward, although undoubtedly the corresponding result must be very complicated 
due to the spin dependence of the low energy interacting Hamiltonian. The spin dependent 
interactions on the light-front are essentially related to the chiral symmetry breaking. These 
spin dependent interactions in HQET are suppressed in the leading order approximation 
because they are l/mq corrections and can be treated perturbatively with respect to the 
heavy hadron states. This is why for heavy quark systems the chiral symmetry breaking 
can be treated separately from the confinement. 

In fact, the model-based theoretical investigations on heavy quarkonia lasted for one 
and half decades is recently replacing by first-principles exploration on QCD. The lattice 
QCD simulation may give an acceptable description for heavy quarkonium spectroscopy with 
manageable control over all the systematic errors The development of nonrelativistic 
QCD provides a general factorization formula to quarkonium annihilation and production 
processes so that a rigorous QCD analysis may become possible ^ . Meanwhile, in the past 
five years considerable progress has been made for heavy hadrons with one heavy quark, 
due mainly to the discovery of the so-called heavy quark symmetry (HQS) [|l^ and the 
development of the heavy quark effective theory (HQET) from QCD. The HQS and 
HQET have in certain contents put the description of heavy hadron physics on a QCD- 
related and model-independent basis. Yet, a truly first-principles QCD understanding of 
heavy hadrons is still lacking since no good nonperturbative QCD approach is available for 
a direct computation of heavy hadron wavef unctions. On the other hand, in the last decade, 
the investigations of the light-front field theory on nonperturbative bound state problems 



have made some progress |jTT[ but no real hadronic problem has been solved from which. 
Starting with heavy hadrons may provide a possible explicit solution of hadronic bound 
states in light-front QCD. 



Simply speaking, the approach to achieve the QCD description of hadronic bound states 
that we shall study in this paper can be summarized as follows: Applying the similarity 
renormalization group approach to light- front QCD, we can obtain an low energy QCD 
Hamiltonian which is an expansion in terms of the QCD coupling constant. Then we at- 
tempt to solve from this low energy QCD Hamiltonian the strong interacting bound states 
as a weak-coupling problem. The weak-coupling treatment contains the following steps: (i) 
Compute from the similarity renormalization group scheme the low energy Hamiltonian Hx 
at the low-energy cutoff A up to the second order in coupling constant. Then separate the 
Hamiltonian into a nonperturbative part, Hxq which contains not only the free Hamiltoni- 
ans of quarks and gluons but also the dominant two-body interactions, and the remaining 
part plus the higher order contributions generated in the similarity renormalization as a 
perturbative term, Hxi- (ii) Introduce a constituent picture which is an important step in 
the realization of the weak-coupling treatment of nonperturbative QCD. The constituent 
quarks and gluons have masses of a few hundreds MeV, and these masses are functions of 
the cutoff A that must vanish when the effective theory goes back to the full QCD theory. 

(iii) Solve hadronic bound states with Hxo nonperturbatively in the constituent picture and 
determine the cutoff dependence of the constituent masses and the coupling constant. The 
coupling constant g now becomes an effective one, gx- In the nonperturbative study of Hxq, 
if we could show that with a suitable choice of the low energy cutoff A, the effective couphng 
constant gx is arbitrarily small, then a weak-coupling treatment could be applied to the low 
energy QCD Hx such that the corrections from Hxi can really be computed perturbatively. 

(iv) There should be a limit gx ^ Qs, where gs is the fixed physical coupling constant mea- 
sured at the hadronic mass scale, such that all the constituent quarks and gluons become 
current ones again. Then the effective low energy theory returns back to the full QCD the- 
ory. If everything listed above works well, we arrive at a true weak-coupling QCD theory of 
the strong interaction for hadrons. 

In the following, I begin with a general bare interacting Hamiltonian with a high energy 
cutoff A that removes the usual ultraviolet (UV) divergences. Then using the light-front 
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similarity renormalization group, we construct a Hamiltonian under the low energy cutoff 
A. The low energy cutoff is introduced via a smearing function in the similarity renormal- 
ization group that effectively integrates over all the modes above the cutoff A. The choice 
of the smearing function in this paper is much simpler in comparison to the original setup 
H. Applying this general formula to the light-front heavy quark effective theory (light-front 
HQET) we developed recently , a low energy confining QCD Hamiltonian can be explic- 
itly obtained for heavy hadron systems. Consequently, a clear light-front picture of quark 
confinement emerges. Furthermore, based on the idea of weak-coupling approach, we use 
this confining Hamiltonian to study the light-front heavy hadron bound states nonpertur- 
batively, and we can then provide an explicit description of the weak-coupling treatment in 
the light-front HQET for heavy hadrons. 

The present paper is organized as follows. In Section 2, we present a general procedure 
of constructing a low energy light-front QCD Hamiltonian in the similarity renormalization 
group scheme and discuss the possible existence of confining interaction in such a derivation. 
In Section 3, applying the general procedure to the light-front HQET of QCD, we further 
derive the heavy quark confining Hamiltonian at the low-energy scale. In Section 4, the light- 
front heavy hadronic bound state equations in the constituent picture are developed within 
the scheme of the weak-coupling treatment. A light-front picture of quark confinement 
for heavy hadrons is illustrated in Section 5. In Section 6, the heavy hadron bound state 
equation is solved for quarkonia with a Gaussian-type light-front wavefunction ansatz, from 
which the scale dependence of effective coupling constant is determined as a solution of the 
similarity renormalization group equation on the quarkonium binding energy in Section 7. 
The low energy running coupling constant is also qualitatively obtained. In Section 8, a 
connection between the low energy effective theory to the full QCD theory is explored and 
the consistency is provided. Finally, a summary is presented in Section 9. 
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II. LOW ENERGY QCD HAMILTONIAN IN SIMILARITY 
RENORMALIZATION GROUP SCHEME 



We begin with the general formulation of the similarity renormalization group approach 
to construct a low energy QCD Hamiltonian. In general, for a given bare Hamiltonian, 
= Hq +Hf , where Hq is a bare free Hamiltonian and Ei is assumed to be its eigenvalue, 
the similarity renormalization group approach leads to the following Hamiltonian at the low 
energy cutoff A (for a detailed derivation, see Ref. |T|): 

+ [[Hf\", Hfyix] , Hfx'r] + [Hfy, [Hfx", Hfx„j] ] ] + ... 
V " ^ 'R' " 'T' / 

" ^ 'R ' « 'R 

= Hf^ + H^x^ + iff + . . . , (2.1) 



where H^,^ = fxijHP^ (we use the notation A, = {i\A\j)), Rfx^.j = -ePe- [ixf^ij) ^up 
and 



X'ij 

'R 



poo 

~fxij / dX'Xx'ij, (2.2) 

J X 

1 / d \ r°° 1 

= ~ ^, _ ^, V^^'^^i) Jx '^'^'^A'ij + ^ _ ^ (1 - fxij)Xxij. (2.3) 

The function fxij = f{xxij) is a smearing function in the similarity renormalization group, 
and Xxij = ^]~^^x ■ "^^^ smearing function is introduced to force the Hamiltonian Hx 
becoming a band diagonal form in energy space. This requires the following properties for 
fxij'. when x < 1/3, / = 1; when x > 2/3, / = 0; and / may be a smooth function from 1 
to for 1/3 < X < 2/3. Thus, through the similarity renormalization group, we eliminate 
the interactions between the states well-separated in energy and generate the Hamiltonian 
of eg. (|2.1| ). The expansion of eq.( |2.1| ) in terms of the interaction coupling constant brings 
in order by order the full theory corrections to this band diagonal low energy Hamiltonian. 

Explicitly, the bare Hamiltonian input in the above formulation can be obtained from 
the canonical Lagrangian with a high energy cutoff that removes the usual UV divergences. 



8 



For light-front QCD dynamics, the bare Hamiltonian in our consideration is the canonical 
light-front QCD Hamiltonian that can be either obtained from the canonical procedure in the 
light-front gauge p!2| , p!3| or generated from the light-front power counting rules Instead 



of the cutoff on the field operators which is introduced in ref . |1| , we shall use in this paper 
a vertex cutoff to every vertex in the bare Hamiltonian: 

e{K'/P+-\pr-p-\)^ (2.4) 

where p~ and are the initial and final state light-front energies respectively between the 
vertex, A is the UV cutoff parameter, and the total light-front longitudinal momentum 
of the system we are interested in. Eq.( |2.4] ) is also called the local cutoff in light- front 
perturbative QCD ||19[. All the A-dependences in the final bare Hamiltonian are removed 



by the counterterms so that the bare Hamiltonian used in eq.( p.lD has already been 
renormalized as A — oo. The use of eq. (|2.4|) largely simplifies the analysis on the cutoff 
scheme in ref. 

Meanwhile, in similarity renormalization group calculation, we should also give an ex- 
plicit form of the smearing function fxij- One of the simplest smearing functions that satisfies 
the requirements of the similarity renormalization group scheme is a theta-function: 

h^J = 0{\-^^^3)■ (2.5) 

However, using the definition of x\ij, we can further replace the above smearing function by 
the following form on the light-front: 

/a.,=^^(;^-|AP,7|), (2.6) 

where A is a low energy cutoff, and AP^y = — P^^ is the light-front free energy difference 
between the initial and final states of the physical processes. The light-front free energies 
of the initial and final states are defined as sums over the light-front free energies of the 
constituents in the states. The smearing function eq.( ^.6D satisfies the requirements for the 
similarity renormalization group approach although it is not a smooth function. 



Throughout this paper, we shall always use the definition of eq.( |2.6| ). Thus, the Hamil- 
tonian (pTT|) can be reduced to 



The front factor (the theta-function) in the above equation indicates that the Hamiltonian 
Hx describes the low energy interactions (with respect to the cutoff A). Therefore, A should 
be a value of the hadronic mass scale. Eq. (p.7|) shows that the low energy Hamiltonian is 
apparently a modified Hamiltonian perturbative expansion. The function g\ij in eq. (|2.7|) is 

9M,k = j^^^^^ d{\'/P+)h>.u^^^jTJ^^fy^k = ei\APr^\ - \yP^)9{\AP-,\ - \APr^\). (2.8) 

The theta function in eq.( p.8|) guarantees that the singularity coming from the small energy 
denominators in the usual Hamiltonian perturbation theory does not occur in the above 
formulation. 

Up to this point, we have introduced two cutoffs, A and A. The cutoff A is used to 
remove the UV divergences as in the usual perturbation theory. While the low energy 
cutoff A is introduced to separate the high and low energy state interactions, it is indeed a 
scale parameter described by the light-front similarity renormalization group. The similarity 
renormalization group transformations integrate out all physical degrees of freedom above 
the cutoff A and generate the low energy Hamiltonian Formally, when A ^ oo, = 1, 
and eqs. (|2.2| ) and ( p.3|) vanish so that Hx — > H^. In practice, we shall take A ~ 1 GeV. 
Of course, the final physical observables must be A independent. As a consequence of this 
condition, the coupling constant and the constituent masses in the Hamiltonian Hx become 
functions of this low energy cutoff A. It is these cutoff (or scale) dependences that link the 
effective theory to the full theory in the limit g —>■ gs- 

Next, to calculate explicitly the low energy effective QCD Hamiltonian, we begin with 
the canonical QCD theory. It is convenient to use the light-front two-component formulation 
of canonical QCD in the light-front gauge A'^ = 0. The bare light-front QCD Hamiltonian 
in such a two-component formalism is given by [jl3 
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= Jjx+d^x^(Ho + Hi). 
Here means that the local cutoff, eq.( p.4| ), has been imposed, and 



(2.9) 



T-ti = Hqqq + 'H„„„ + 7^„„„„ + + H„„„„ + countcrterms 



(2.10) 

f2.in 



with 



'H,,, = 9^ [-2 (^) {d-A^) + a. (ip) (a ■ 9^ + m) 



n 



... - gn \ d^AiAlAi + {d^AD [ — ) {Aid^Ai) 



1 



+ 2 (^) (/'^^'=A^,9+A^) (i^ ) (^tT'^^) 



^-- = 2/{(^) (^'^''O (^) (m)} 



.... 



+2 ( i^) {J^,d^A\) (i, ) 



(2.12) 
(2.13) 

(2.14) 
(2.15) 

(2.16) 



where A^ = T"'Aa± is the transverse component of the gauge field, T° is the generator of 
SU(3) color group, and ^ is the two-component form of the light-front quark field: 



^_ = A ^ 







(2.17) 



The notation cr is defined by: = cr^, cr^ 
use of the light-front 7-representation: 



(•Jt) [a'{td^ + gA')+tm]^ 
(the Pauli matrices). This comes from the 
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, 7^ = 
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, 7^ = 



















f2.18) 
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which is different from the one we used in ref. |T^. The reason of using this new 7- 
representation is that it leads to not only a realization of the two-component form for 
the light-front fermion field, but also a correct correspondence of the fermion spin opera- 
tor 5** ~ (T* on the light-front. Counterterms are added to eq.( |2.11| ) in order to remove 
all the A-dependence. Thus the coupling constant g in eqs.( |2.12| - [^.16D is a perturbatively 
renormalized coupling constant. 

Upon the calculation to the second order in the coupling constant, we obtain an effective 
Hamiltonian (|2.7|) in qq sector (with the initial and final states, \i) = b'^ {pi, Xi)d^ {p2, \2)\0) 
and \j) = &^(p3, X3)d^{p4,, A4)|0), respectively, where pi and Aj denote the respective momen- 
tum and helicity of a quark on the light-front): 



H,,,=9{\yP^-\APr\)l{J\■.H^■.\^) 



1 

P+ 



9 

X^Cf-;^ Ine + mass counterterms 



-g\T'^){T 



27r2 

,qi ^-Ps 



im , 



q^ 



[pi 



-i'O" ■ pi± + ^""^ 

a - 



,ql 



a ■p2± +im ^■ 
-a 



a ■ 



[pt] . 



(2.19) 



9^ [Pi \ [pt] 

Here P"*" = pt + pt = P3 + ? ^^ij = Pi P2 ~ Ps ^ Pa ^ '■ : represents a normal 
ordering, where the instantaneous interaction contribution to the quark self-energy has been 
included in the self-energy calculation which is given by the mass counterterms and the 
logarithmic divergence in (g]Tg), the color factor Cf = {T'^T'') = (iV^ - 1)/2A^, = 3 the 
total numbers of colors, and e is an infrared longitudinal momentum cutoff. Since Ine is 
an infrared divergence, it cannot be removed by mass counterterms. In gauge symmetry, 
this divergence must be canceled in the physical sector (and this is true as we will see 
later). The last term in (|2.19| ) is the one-gluon exchange contribution to the low energy 
Hamiltonian. The momentum q is carried by the exchange gluon: q'^ = pt ~ Pt = pt ~Pti 
q± = pi± — p3± = p4± — p2±, Xa, denotes a helicity eigenstate. The factor Frij arises from 
the similarity renormalization group scheme: 



rij 



0(aVp+ - i^r -P3" - q-[)9{AVP^ - [PI-P2 - q-\) 
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^(bi -Ps -Q I - p+ )^(bi -P3 -Q \-\P4 -P2 -Q I) 



Pi -P3 -Q' 

A2 



^ djlPi -P2 -Q I - -p+Wi\Pi -P2 -Q \-\Pi -P3 -Q I) 



P4 -P2 

{pi < — ^ P3 , P2 < — ^ P4 , g — ^ -q)}- (2.20) 



All light-front energies in eq.( p.20|) are on mass-shell: = ^'^^+^' and q ^ 



If we continue to evaluate all the terms in the expansion of eq.( p.l| ), the resulting Hamil- 
tonian is the exact QCD Hamiltonian. In practice, we should only consider the leading and 
the next-to- leading terms, i.e., eq. (p.l9| ), as a starting effective Hamiltonian. The basic idea 
to realize a weak-coupling treatment of QCD for hadrons is whether we can solve hadron 
states from this effective Hamiltonian ( |2.19| ) with an arbitrary small coupling constant g 
such that the higher order corrections in ( ^.11 ) can be handled perturbatively. From the suc- 
cess of the constituent quark model, we understand that a necessity for such a realization is 
the existence of a confining interaction in (|2.19|) . 

Naively, we know that any weak-coupling Hamiltonian derived from QCD will have 
only Coulomb-like interactions, and confinement can only be exhibited in a strong-coupling 
theory. However, Perry has recently found that upon to the second order calculation of 
the low energy Hamiltonian, a logarithmic confining potential has already occurred [Q. 
Explicitly, the two-body quark- ant iquark interaction from the first term of eq. ( |2.19|) is the 
instantaneous gluon exchange interaction 'Hqqqq which has the form in the momentum space: 

The confinement must be associated with the interaction where g"*" — > 0. This is because only 
these particles with zero longitudinal momentum can occupy the light-front vacuum state 
Thus, we need only to analyze the feature of the effective Hamiltonian when g+ 0. 
For g+ ^ 0, the dominant second order contribution from one transverse gluon exchange 
interaction in eq.( p.l9| ) is given by 



q\ q 



e{q^-\yP+). (2.22) 
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In the usual perturbative calculation, no such a ^-function is attached in the above equation. 
Thus these two dominant contributions from the instantaneous and one-gluon exchange 
interactions are exactly cancelled when q'^ —>■ 0. Only Coulomb-type interaction (the terms 
~ 0{-^)) remains. However, in the light-front similarity renormalization group scheme, 
the one-gluon exchange contribution to eq. p.22| ) only contains these gluons with energy 



being greater than the energy cutoff A^/P^. As a result, the instantaneous gluon exchange 
term l/(g"^)^ remains uncancelled if the gluon energy, qj_/q'^, is less than A^/P^. The 
remaining uncancelled instantaneous interaction contains an infrared divergence and a finite 
part contribution (for a detailed derivation, see section V). The divergence part is cancelled 
precisely for physical states by the same divergence in the quark self-energy correction [see 
( P.19| )]. The remaining finite part corresponds to a logarithmic confining potential: 

bxln\x~\ + Ca In ) . 

The above result is first obtained by Perry with the use of the concept of coupling coherence 
and a slightly different renormalization scheme |^ (also see a oversimple derivation given by 
Wilson ) . Here the derivation is purely based on the light-front similarity renormalization 
scheme 

One may argue that the existence of such a confining potential in Hx may only be an 
artificial effect designed in the renormalization scheme we used. If we included the inter- 
action with the exchange gluon energy below the cutoff, then the instantaneous interaction 
would be completely cancelled, and no such confining potential should exist, as expected 
in the usual perturbation computation. Wilson has pointed out that the set up of the new 
renormalization scheme is motivated by the idea that the gluon mass must be nonzero in 
the low energy domain (a constituent picture), which could be regarded as an effect of the 
nontrivial low energy gluon interactions. The cutoff A is of the same order as the constituent 
gluon mass. Thus the gluon energy cannot run down below the cutoff A. The existence of 
the confining potential is a result of the low-energy gluon interactions. In contrast, the 
photon mass in QED is zero at any scale. The instantaneous photon exchange interaction 
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is always cancelled by the corresponding one transverse photon exchange interaction, and 
only the Coulomb interaction is left. Thus, the above confining potential can only exist in 
QCD [|T^. More detailed discussions will be seen in the subsequent sections. 

Yet, even up to the second order, the effective QCD Hamiltonian H\ is already very 
complicated. In order to to examine the above ideas of confining mechanism and to develop 
explicitly a weak-coupling treatment approach of nonperturbative QCD to hadronic bound 
states, in the next section we shall utilize the above formulation to heavy quark systems. We 
find that the low-energy QCD Hamiltonian H\ for heavy quarks can be largely simplified 
and an analytic form consisting of the confining potential and Coulomb potential emerges. 

III. LOW-ENERGY HEAVY QUARK CONFINING HAMILTONIAN 

In the past few years, QCD has been made a numerous progresses in understanding the 
heavy hadron structure, due mainly to the discovery of heavy quark symmetry by Isgur and 
Wise [|r^, and the development of the heavy quark effective theory by Georgi et al. Very 
recently, we have reformulated the heavy quark effective theory from QCD on the light-front 
PJ^, which may provide a convenient basis for the further study of the nonperturbative 
structures of heavy hadron. Now we use the light-front HQET to derive the low energy 
heavy quark effective Hamiltonian in the similarity renormalization group scheme. 

A. Light-front heavy quark effective theory 

The light-front heavy quark effective Lagrangian derived from QCD lagrangian C = 
Q{ip — mQ)Q as a l/mg expansion is given in refs. 

9 OO -1 

C = —Ql+{iv ■ D)Q,+ - J2 ( -YoJ^+Uia ■ D){-iD+Y-\ia ■ (3.1) 

where 2^+ is the light-front dynamical component of the heavy quark field after the phase 
redefinition: 

Q{x) = e-'^<^^-%Q,+ {x) + Q,^{x)), (3.2) 
15 



v'^ the four velocity of the heavy hadrons, = Mhv^ with f ^ = 1 and Mh being the heavy 
hadron mass, a- D = a± - — -^{oi±_ ■ w_l + P)D~^, and is the usual covariant derivative. 
The corresponding light-front heavy quark bare Hamiltonian density is given by 

+ E ( TTQl+k"^ ■ D){-zD^r-\za ■ . (3.3) 

In the large mg limit, only the leading (spin and mass independent) Hamiltonian is 
remained. In other words, the phase redefinition ( p.2| ) removes the dominant piece of the 
space-time dependence of the heavy quark. The remaining dependence is only due to the 
residual momentum of the heavy quark in the heavy hadrons. The l/rriQ terms {n > 1) in 
( p.3| ) can be regarded as perturbative corrections to the leading order operators and states. 
Therefore, the heavy quark mass rriQ is indeed a factorization scale for separating heavy 
quark short and long distance dynamics. To determine confining interactions between two 
heavy quarks or a heavy quark with a light quark, only the leading heavy quark Hamiltonian 
plays an essential role. We choose the light-front gauge ^4+ = 0, the leading-order bare QCD 
Hamiltonian density (corresponding to the limit of niQ oo) is given from (|3.3|) : 



'Hid '■ 



1 



QU{v^ (^)d^-A^\-v^-A^}Q, 



T^O + Hqqg + Ti.qqqq. (3.4) 



Note that besides the leading term in eq.(|3.3|), the above bare Hamiltonian has already 



also included the relevant terms from the gauge field part, — |Tr(F^,^F^'^), of the QCD 
Lagrangian. These terms come from the elimination of the unphysical gauge degrees of 
freedom, the longitudinal component A~, in the light-front gauge (see a detailed derivation 
in refs. |ll2|JT3[| ). Eg ( p. 41) has obviously the spin and flavour heavy quark symmetry, or simply 
the heavy quark symmetry. 
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In momentum space, the free part of the bare hght-front effective heavy quark Hamilto- 
nian can be simply expressed by 

^0 = E / [d'k]^{2v^ ■ - v-k+){bl{k, X)K{k, A) + dl{k, X)d,{k, A)} 

^ J V 

= E / [d^k]k-{bl{k, X)K{k, A) + dl{k, X)d,{k, A)}, (3.5) 

where k is the residual momentum of heavy quarks, p'*' = rnqv'^ + k^, and A its helicity. We 
have introduced the notation for the space components of hght-front momentum {p^ iPi) = p 
so that [d^p] = '^2(2'tp'^ • '^^^ operator h\{k,X) [(iJ,(A;, A)] creates a heavy quark [antiquark] 
with velocity f , residual momentum k and helicity A, 

{K{k, A), hl,{k', A')} = {d.{k, A), dl,{k', A')} = 2{2T:f5,,,5xx'5\k - k'), (3.6) 

where 6^{k - k') = 6{k+ - k'+)6\k^_ - k'^). 

Eq. (|3.5| ) means that after the redefinition of the heavy quark field, the heavy quark in 
the light-front HQET carries the effective free light-front energy 

k- = ^{2v^-k^-v-k+). (3.7) 

The meaning of this result becomes more transparent if we expand the light-front energy 
dispersion relation of the heavy quarks as an inverse power of mg: 



p~= ^ p+ = ^Q^~ + :^(^^-L ■ ~ v~k'^) + 0{l/mQ) 



rnqv + k . (3.8) 



We see that the mass part {mQV~) has been removed by the phase redefinition of ( p. 21) . 
Thus, eq. (|3^) describes effectively the "lighten" heavy quark dynamics with respect to its 
residual momentum. In other words, the heavy quarks in the effective theory have the same 
energy scaling behavior as the light ones. 

The above leading Hamiltonian (or Lagrangian) is the basis of the QCD-based description 
for heavy hadrons containing a single heavy quark, such as B and D mesons. However, as 
recently pointed out by Mannel et al. [|l^,|l3l the purely heavy quark leading Lagrangian 
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may be not appropriate for the description of heavy quarkonia states. This is because 
the anomalous dimension of the QCD radiative correction to the QQ currents contains an 
infrared singularity in the limit of two heavy constituents having equal velocity. Such an 
infrared singularity is a long distance effect and should be absorbed into quarkonium states. 
To avoid this problem, they argued that one may incorporate the effective Hamiltonian with 
at least the first order kinetic energy term into the leading Hamiltonian fl^. The kinetic 
energy in light-front HQET is given by 

n,,n = — ^—QlJdl - + ^-oAq.^. (3.9) 

As a consequence, in the heavy mass limit, quarkonia have spin symmetry but no flavour 
symmetry. In momentum space, 

Hun = E / [d'k] (kl-2v^ ■ kX + 

rnQV^ ^ J ^ f+ f + ^ 

^[hl{k,\)K{k,\) + dl{k,\)d^{k,\)]. (3.10) 

The kinetic energy of ( |3.10| ) can be simply obtained by expanding ( p.8| ) up to order l/m^. We 
will discuss later the effect of this kinetic energy in the determination of heavy quarkonium 
bound states. 



B. Low-energy effective Hamiltonian for heavy quarkonia 

Within light-front HQET, we can follow the procedure described in the previous section 
to find the effective low energy QCD Hamiltonian for QQ systems. The bare Hamiltonian 
for QQ systems is given by ( p.4|) plus the leading kinetic Hamiltonian (|3.9| ) for both heavy 
quark and antiquark, where two heavy constituents have the same velocity carried by the 
heavy quarkonia. The kinetic energy which is of order Kqcd/^Qi is at most the same order 
as the Coulomb interaction. It may affect on the quarkonium bound states but not on the 
derivation of the long distance quark interactions. In fact, the success of the potential- 
model description indicates that the scalar interactions between the two heavy constituents 
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are flavour-independent. In other words, the confining interaction in QQ states should be 



independent of the kinetic energy (|3.9| ). Therefore, we may treat the kinetic energy as the 
same as the instantaneous QQ interaction [the last term in eq.( p.4|) ]. In the derivation of the 
low energy heavy quark Hamiltonian, the free Hamiltonian used in similarity renormalization 
group scheme is then simply given by eq.( |3.5| ). 

With the above consideration, it is easy to find that the leading order contribution to 
the low energy effective Hamiltonian is the low energy part of the heavy quark effective bare 
Hamiltonian, 
A2 



H 



(0) 



p+ 



|Ai^-|)(j| j dx-<fx^\^-^Ql+{v-d+ - 2v^ ■ d^)Q,, 



1 



19, - A, 



+ 



- + ^9+2 \\^), (3.11) 



plus all other l/mg terms in ( |3.3| ) as well as the light quark and gluon full QCD Hamiltonian 
that has not been included in the above equation [see (|2.9|) ]. Here the initial and final states 
are defined by \i) = bl{ki, Xi)dl{k2, \2)\0) and \j) = &J,(/c3, A3)(iJ^(A;4, A4)|0), respectively, 
= Pi + P2 = {^Q + ^q)'V^ + kf + k2 = {tuq + mQ)v~^ + k^ + 

(2) (2) (2) 

The next-to-leading order contribution contains two different parts, ^ = H^-^^ + -ff]^2 ? 

(2) 

where is the self-energy correction, 



[£k]' 



' e{ki - k+) /{ki - k)i 

kt -k+ \ kt - k+ ~ ^ 



X 



|fc--A:--(fci-fc)-|-^) 
ki — k~ — {ki — k)~ 



k--{k^-k)-\) 



+ 



e{k+ - k+) f {k2 ~ k)i V 



■J N 2 



— k+ 

rvn fx, 



k+ 



k--{k2-kr\) 



1^2 ^ 



{k2 - k)- 



p+' 



{k2 - k)- 



e{— - \AP,r^m2nYYS%k, - k,)6\k2 - k,)6,,,,,6x,x. 
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^ J 2(2vr)3 ix-xy ^'' - '^^ - ""'^ 



X (3.12) 



(2) — 

and H^2 is the QQ interaction, 



X 



[ g(|fci - ^3 - q-\ - ^)e{\k^ - ^3 - g~| - - fca - g' 
1 K - K - q- 



+ 



0{\ki -q-\- ^)e{\k^ -k2 -q~\-\h " ^3' " r I) 

^4 ^ ^2 "9 



= ^^(-^ - I AP,7|)2(27r)353(fc^ + ^2 - ^3 - ^4)5Aa.<^A.A4i^(a: - x', - n'^, Mh) 

(P+)2 _ a;/)2 _ ^/^)2 + _ x')2mI ' ^ " ^ 

where /c~ is given by eg. (|3^) , q^ = kf — k^ = kf — k2 , q± = — ^3^ = ki± — k2±, and 
q^ = q\/q'^ , we have also introduced the longitudinal residual momentum fractions and the 
relative transverse residual momenta, 

X = kf/P^ , K± = ki± — xP±, 

x' = k^/P+ , k'^ = k3i_ - x'P±, (3.14) 
and defined the function F, 

F{x, k, M) = e{A{x, k, M) - X'^)e{A^ - A{x, k, M)), (3.15) 

with 

k^ 

A(x,k,M) = — + \x\M^. (3.16) 

|x| 

Since < Pi = mqv'^ + k^ < P+ = MfjV^, in the heavy quark mass limit, we have 
Mh — i> 2mQ so that —rnqv^ ^ k^, k^ < mqv^ . Hence, the range of x and x' is given by 

-l<x, x'<i. (3.17) 
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Eqs. (|3.11| ), ( p.l2|) and ( |3.13| ) consist of the effective Hamiltonian for quarkonia up to the 



second order in the similarity renormahzation group scheme. 

Apparently, the above effective Hamiltonian is not a low energy Hamiltonian because 
P+ = Mhv^ which is of order a few GeV. However, from eq. ( |3.5| ) and eq.( |3.8| ) we see that in 
the bare heavy quark Hamiltonian the mass term rnqv" has been integrated out. To address 
the correct energy scale of the low energy heavy hadron dynamics, we should introduce the 
residual center mass momentum of the heavy quarkonia and the residual heavy hadron 
mass A, 

K>' = Av'' , A = MH-mQ-mQ, (3.18) 
where v'^ is still the four-velocity of hadrons. It follows that 

K-^ = kt + kt = kt + kt , K^ = ki^ + k2± = ks^ + ki±. (3.19) 



With the residual heavy hadron momentum K'^ considered, eqs. (|3.12|) and ( p.l3 ) become 
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H^x^ = - \AK.^\m2n)rS'ik, - k,)5\k, - k,)5,,,,5,,,, 

^ J 2(2vr)3 (y-yy ^y " '^^ " 

X _ (3.20) 

hS = - \AKT-^\)2{27:f5\h + ^-1,- h)6x,x,6x,x,Fiy - y\ - k'^, a') 

V(T")(T") 1 {n^ - n'^f 

{K+Y {y - y'f (;,^ _ t^'^)2 + _ y^yx^ ' 

where we have also introduced the corresponding residual relative momenta: 

y = kf/K^ , K_L = ki± - yK±, 

y' = kt/K+ , k'^ = ks^ - y'K^. (3.22) 

The range of the residual longitudinal momentum fractions y and y' are given by 

— oo < y = < oo , — CO < y = < oo. (3.23) 

Now all the quantities appearing in the effective Hamiltonian have the low energy scale of a 
few MeV. 
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C. Reexpression of the low energy Hamiltonian in the weak-couphng treatment 

scheme 



As we have mentioned in the Introduction, the first step to follow the idea of the weak- 
coupling approach is to construct the low energy Hamiltonian H\ up to the second order 
and then separate it into Hxq and Hxj, 

Hx = Hxo + Hxi. (3.24) 



In eq.( |3.24|) , if^o contains the free Hamiltonian plus the dominant two-body interactions 



which conserve the particle number, and Hxi is the remaining interaction Hamiltonian which 
describes the emission and reabsorption processes plus all the higher order terms in the ex- 
pansion of eq.(p.l|). Once H\ is derived, we can reexpress it as eg. ( \3. 24\ ) such that H\q is 



set up to nonperturbatively determine the hadronic bound states, and H\i should be treated 
perturbatively. This separation is a basic step to realize a weak-coupling treatment of non- 
perturbative QCD [|l|. Thus, besides the free quark and gluon Hamiltonian with constituent 
masses, H\o also contains the instantaneous interaction and all the second order contribu- 
tions generated by integrating over all the modes above the low energy cutoff A, namely, 
eqs.(F^ and ([Oil ). Explicitly, 



H,o^, = 9{\yK^ - |Air^|){iJQQ^_ + Vqq,}, (3.25) 

where 

HQQfree= [2(2^)'? 5' (kl " h)S'ik2 ' h)Sx,xA.X. 

X lk^ + k^ + -^ + -^-2^Cf^\ne], (3.26) 

1 ^ 2 2mo 2m7T Att^ ^ K+ ' ^ ^ 



3.3n T T . -Ag\T''){T'') \ 1 

I \2 



{K+Y \{y-y'y 



""^^ -,eiAiy-y',K^-K'^,A)-X') 



(y - y'y (kx - /tl)2 + iy- y'YA 
2{2Txf5\ki + h-h- h)Sx,x.A2M 
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-r-2 

+ ^^^(^(Z/ - ^± - '^l, A) - A^) . (3.27) 

{k± - K'_^y + {y- y'yA ) 

In (|]26]), k;- is given by (^), kf = ^{k\- j^2v^ ■ k±k+ + ^^k+^). In the above results, we 

have already let UV cutoff parameter A ^ oo and the associated divergence has been put in 

the mass correction. The Icinetic energy ( p.lO|) now is included in the above nonperturbative 



part of the effective Hamiltonian. 

The free energy part -f^ggj^ee eq.( |3.26| ) has also included the self-energy correction 



which is the instantaneous interaction contribution, a normal ordering term of the instan- 
taneous interaction in ( |3.11D , plus the one-loop contribution ( p. 201 ). The result is 



S = Vc 



dyid^Ki^ \e{y - yi) 



2(27r)3 \{y-y,Y 



^F{y - yi,K±- Ki±,A) 



2 f dyid^Ki±e{y -yi) ^ 
^ J 2(27r)3 iy-y,y ^ ~ A{y - y^, - A)) + 6mQ 

= -£^X'Cf\ne + 6ml, (3.28) 

where the mass correction SruQ = ^C/A^ln^ which has been renormalized away in 
eg. (13.261). By removing away this mass correction, we should assign the corresponding 



constituent quark mass in Hxq being A-dependent. But, the heavy quark mass is much 
larger than the low energy scale. Its dependence on A should be very weak and could be 
neglected. 

The QQ interaction Vqq of eq. ( p.27|) contains the one gluon exchange interaction 



eq. (|3.21|) plus the instantaneous interaction [the last term in ( p.ll| )]. It clearly shows that 



without the low energy cutoff (A = 0), the instantaneous interaction is completely cancelled 
by the same contribution from the one transverse gluon exchange and the remaining one 
gluon exchange interaction is a Coulomb interaction, like in QED. Now, with the low en- 
ergy cutoff, the one gluon exchange contribution only contains these gluons with the energy 
greater than the cutoff A. Thus, the resulting QQ interaction has two terms: The first term 
is the result of the noncancellation between the instantaneous interaction and one trans- 
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verse gluon exchange interaction, which corresponds to a confining potential. The second 
term, the rest of one transverse gluon exchange interaction, is the Coulomb interaction on 
the light-front. The detailed confinement mechanism on the light-front will be discussed 
in section V. With the kinetic energy incorporated, we see that the above effective QCD 
Hamiltonian which will be used to determine the heavy quarkonium bound states only has 
the spin symmetry but no flavour symmetry. 

Before ending this section, we may compare the present formulation for heavy quarkonia 
with the nonrelativistic QCD formulation developed by Lepage et al. 

In the nonrelativistic QCD formulation, heavy quarkonia are described by an effective 
fleld theory of QCD in the nonrelativistic limit plus a systematic computations of the rel- 
ativistic corrections (in terms of momentum scales (Mf)^ and {Mv^)"^) and QCD short 
distance corrections (in terms of the scale M^). The nonperturbative QCD scale Aqq,j^ is 
implicatively included in this formulation. Heavy quarkonium annihilation and production 
processes can then be factorized with respect of the above different scales. The domi- 
nant contributions in quarkonium processes, namely, the nonperturbative QCD dynamics of 
quarkonium bound states may be computed in lattice simulations. 

Our formulation is based on the factorization scale tuq which naturally separates the 
QCD short distance and long distance dynamics. The long distance dynamics is described 
by the residual momentum which is now controlled by A}^cd/ effective Hamiltonian 
Hxo- The resulting effective Hamiltonian derived from QCD by the similarity renormaliza- 
tion group approach contains explicitly the conflning and Coulomb interactions which have 
encompassed the necessary long distance effects for heavy quarkonia. The quarkonium bound 
states can then be directly solved in the corresponding light-front bound state equation (as 
we shall see later). The short distance dynamics can be systematically computed in the 
ordinary perturbation theory. These include the QCD radiative corrections (controlled by 
ln(mQ/A), where A is an UV cutoff), and the l/mg corrections (controlled by A/mg). Our 
formulation is fully relativistic. It is a more complete QCD formulation for heavy quarkonia 
in comparison to nonrelativistic QCD [^. It allows to directly compute the nonperturbative 
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QCD dynamics without the help of lattice simulation. Moreover, it is also straightforward 
to extend this formulation to the heavy hadron system which contains a single heavy quark, 
as we shall see the next. 



D. Low-energy effective Hamiltonian for heavy-light quark systems 

The heavy-light quark system (heavy hadrons with one heavy quark) is one of the most 
interesting topics in the current study of heavy hadron physics. We now apply the similarity 
renormalization group approach to such system. 

The bare cutoff Hamiltonian we begin with for heavy-light quark system is the combi- 
nation of the heavy quark effective Hamiltonian (|3.3| ) and the full Hamiltonian for the light 
quark ( p.9|) . Due to the HQET, we may also introduce the residual center mass momentum 
for heavy-light systems, 

= Av+ = p+ + fc+ = p+ + = Av± = pi± + = p2± + fcsx, (3.29) 

where A = — niQ, pi and p2 are the light antiquark momenta in the initial and fi- 
nal states respectively, and ki and k2 the residual momenta of the heavy quarks. The 
initial and final states in Qq sector are denoted by \i) = bl{ki, \i)S{pi, \'i)\0) and 
\j) = 6j(A;2, A2)(i^(p2, A2) |0), respectively. The residual longitudinal momentum fractions 
and the residual relative transverse momenta are defined in the similar way as in quarko- 
nium system, 

y = Ptl , f^± = Pi± - yK±, 

y' = pt/K+ , k'^ = p2± - y'K^, (3.30) 

but the range of the longitudinal momentum fractions y and y' is different: 

0<y = ^—<oo , 0<'i/=^— <oo. (3.31) 

Following the general procedure described in the previous subsection, it is easy to find 
that the nonperturbative part of the low-energy effective Hamiltonian for heavy-light quark 
systems, which is given by 
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H 



XOij 



K+ 



~\^K7^\){HQ7ifree + VQ-,i], 



(3.32) 



where 



Qqfree= [2 (27r) ^'^'^ (^1 " k2)S%Pl - P2)Sx,X,6x.,X^ 



P 



(3.33) 



VQq{y - y\ K± - K,'±) = 2{27vY6\h +Pi~k2- P2)Sx^x,Sx,x, 



with 



A2 A2 

xO{—-\Pl-P2-q-m- 



X 



0{\pi -P2 -q' 



1^2 -ki -Q I) 



K+ 



)Q{\Pi -P2 -Q \-\k2 -ki -q- 



+ 



9{\k2 - fci -q~ 



Pi -P2 -Q' 



a2 



K- 



)9{\k2 -k^ -q \-\pi -P2 -q~ 



2i27ifS%ki+p,-k2-p2)Sx,xA 



^2 ~ ki — q 

-2g\T-)iT-) 



X2X4 ' 



(K^ 



{y - y') 



i\2 



I ^2 



2 I I 

/Sj I ' ^ I ^_L * ^ 



, (/t± - K 

{y-y'y y{y-y') y'{y-y') 
e{B - \^)e{B - A) 

(k_l - k'_l)2 - + m2 



+ 



y-y' /^^2 

^(A- A2)^(A-5) 
(k^ - fi:Y)2 + {y- y'ff\ 



(3.34) 



/ \2 



y-y' 



y 



y' 



A^{^^-^'^r/\y-y'\ + \y-y'\A\ 



(3.35) 
(3.36) 



Here we do not include the heavy quark kinetic energy into Hxq since the dominant kinetic 
energy is given by the constituent hght quark. The heavy quark kinetic energy can be 
treated as a perturbative correction to i^Ao- The heavy quark free energy has been written 
in (PI) by 



k = — ( 2v± ■ k± — V k'^ 



K+ 



2K\ -ki -K'k^ 



(3.37) 
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The low energy heavy-light quark effective Hamiltonian is mg- independent. It obviously 
has the spin and ffavour symmetry, namely, the heavy quark symmetry. Comparing to the 
quarkonium systems, the Qq interactions are much more complicated. But it is not difficult 
to see that the above Hx contains a confining potential. 

Finally, it is also straightforward to extend the above derivation to light-light quark 
systems. The result is just eq.( p.l9|) but in terms of the relative momenta: 



x = pt/P^ , n±=Pi±-xP±, 

x' = P+/P+ , k'^ = pii_ - x'P±. (3.38) 

Here there is no residual center mass momentum for light-light system. The hadron mo- 
mentum P'^ is already of order a low energy scale. We shall not intend to discuss the 
light-light systems in details in this paper. As we have pointed out in the Introduction, 
for the light-light quark systems, besides the confinement, chiral symmetry breaking also 
plays an essential role in the low energy hadronic dynamics. We shall remain the light-light 
quark systems for further investigation when we attempt to address the problem of chiral 
symmetry in light-front QCD [|T^. 

In conclusion, we have obtained in this section the renormalized low energy effective 
QCD Hamiltonian for heavy- heavy and heavy-light quark systems, and extracted the non- 
perturbative part Hxq, eqs. (|3.25|) and (|3.32|) , in the weak-coupling treatment scheme. We 
are now ready to solve heavy hadrons on the light-front. 

IV. LIGHT-FRONT HEAVY HADRON BOUND STATE EQUATIONS 

In this section, based on the low-energy heavy quark effective Hamiltonian derived in the 
previous section, we shall construct light-front bound state equations in the weak-coupling 
treatment scheme. 
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A. General structure of light-front bound state equations in the weak-coupling 

treatment scheme 



In general, a hadronic bound state on the light-front can be expanded in the Fock space 
composed of states with definite number of particles [|T^,^. Formally, it can be expressed 
as follows 

\^{P+,P^,X,)) = Y.j [\{[d^p])2{2T:)H\P-Y.p;)\n,p,\)^n{x,,KuAi). (4.1) 

n,\i i i 

where P~^, P± are its total longitudinal and transverse momenta respectively and its total 
helicity, \n, p, Aj) is a Fock state consisting of n constituents, each of which carries momentum 
pi and hehcity Aj Aj = As); K^j, Aj) the corresponding amplitude which depends 

on the helicities Aj, the longitudinal momentum fractions Xi, and the relative transverse 
momenta K±i: 



Xi 



Pi 
P+ 



K-i± = Pi± - XiP± 



(4.2) 



The eigenstate equation that the wave functions obey on the light-front is obtained from 
the operator Einstein equation P^ = P^P~ — P]_ = M^: 

P2 + M2 



Hlf\P , P±, Xs) 



P+ 



\P^, P±, As), 



(4.3) 



where Hip = P is the light-front Hamiltonian. Explicitly, for a meson wave function, the 
corresponding light-front bound state equation is: 



Xi 



^qq 




$ _ 





{qq\Hint\qq) {qq\Hint\qqg) 
{qqg\Hint\qq) 



qqg 



(4.4) 



where Hint is the interaction part of Hip. 

Obviously, solving eq.( [4.4| ) from QCD with the entire Fock space is impossible. A basic 
motivation of introducing the weak-coupling treatment scheme is to simplify the complexities 
in solving the above equation. The Hamiltonian H\ derived in the previous sections has 
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already decoupled the high and low energy states. Here we only consider the low energy 
states (hadronic bound states). Due to the kinematic feature of boost symmetry on the 
light-front, we can assign a relative small longitudinal light-front momentum to the bound 
states. On the other hand, the light-front infrared divergences force us to introduce a small 
cutoff on the longitudinal light-front momentum to each individual constituent. Thus, the 
hadronic bound states can only consist of the Fock space sectors with a few particles. This 
is a kinematic truncation on eqs. ([4.1| ) and ([4.4|) . Furthermore, the most important point in 
the weak-coupling treatment scheme is the reseparation of the Hamiltonian H\ = Hxq + Hxi. 
As we mentioned before, Hxo which conserves particle number devotes to a nonperturbative 
evaluation to the bound states through eq.( [4.4|) . And Hxi which describes the particle 
emissions and reabsoptions is hopefully a perturbative term in the weak-coupling treatment 
so that we may not consider its contribution to eg. ( |4.4| ). Then, eq.( |4.4| ) becomes diagonal 
in Fock space with respect to the different particle number sectors, 

{qq\Hxo\qq) $ 

{qqg\Hxo\qqg) $ 


Now we see that the bound state equation is manable. 

The second important step in the weak-coupling treatment to the low energy QCD is the 
use of a constituent picture. The success of the constituent quark model suggests that we 
may only consider the valence quark Fock space in determining the hadronic bound states 
from if AO- In this picture, quarks and gluons must have constituent masses. This constituent 
picture can naturally be realized on the light-front However, an essential difference from 
the phenomenological constituent quark model description is that the constituent masses 
introduced here are A dependent. This cutoff dependence of constituent masses (as well 
as the effective coupling constant) is determined by solving the bound states equation and 
fitting the physical quantities with experimental data. This is indeed a renormalization 
condition in nonperturbative QCD. Note that unlike the usual renormalization scheme in 
QED, quarks and gluons in QCD are not physically observable particles so that we can only 



Xi 



(4.5) 
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determine their renormalized masses and coupling constant in hadronic (composite particles) 
sectors. Once the constituent picture is introduced, we can truncate the general expression 
of the light-front bound states to only including the valence quark Fock space. The higher 
Fock space contributions can be calculated as a perturbative correction through Hxi. Thus, 
for mesons, eg. ([4.1D is reduced to the following simple form: 



x$gg(x, K^Ai, A2)|g(pi, A2)g(p2, A2)), (4.6) 



A1A2 ' 



where \q{pi-, \i)q{p2, ^2)) = b^ipi, ^i)d''{p2, ^2)\0), and b^d^ should be regarded as the cre- 
ation operator of the constituent quark and antiquark respectively. Consequently, the con- 
stituent quark and gluon masses and coupling constant g in the effective Hamiltonian 
Hxo become explicit functions of the low energy cutoff A. 

It is worth pointing out that spin is always a troublesome issue in the light-front approach. 
The meson light-front bound state we have constructed is labelled by helicity rather than 
spin. However in practice low-energy hadronic states with definite spins are needed. This 
discrepancy is usually remedied by introducing the so-called Melosh rotation which 



transforms a single particle state from the light-front helicity basis to the ordinary spin 
basis, 

A/(mi + XiMoY + 

where n = (0, 0, 1), and 

= — i + ^ ^ 4.8 
X 1 — X 

With Melosh transformation incorporated, the light-front meson bound state with a definite 
spin can be expressed in the weak-coupling treatment scheme as follows 



A1A2 ' 



|^(P+, P^, J, J,)) = ^ I [d'p,][d%]2{27rf6\P -Pi- P2) 

x$4^^(x, K^, Ai, A2)|g(pi, \i)q{p2, A2)), (4.9) 
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where 

^if{x, K^, Ai, A2) = (j)gq{x, K±)RiX{x, K^), (4.10) 
RiX{x,^±) = E(^il^^(^''^^'"^i)Mi)(^2|i?^(l-a;,-Kx,m2)|s2)(|si^S2|JJ.), (4.11) 

and {^Si^S2\JJz) is the Clebsch- Gordon coefficient. The normahzation condition for the 
state |\I/(f , J, Jz)) is taken to be 

(vl/(P'+,Pl, J', J^)|vl/(P+,P^, J, J,)) = 2{2'nfP+5\P' -P)5rj5j.j^, (4.12) 
which leads to 



dxd K± 



x,K^)\' = l. (4.13) 



2(27r)3 

After the above consideration of the spin property on the hght-front, eq.( [4.4|) , for mesons, 
becomes a hght-front Bethe-Salpeter equation: 



(m2-M2)$4^^(x,«:x,Ai,A2) = (-^X'Cj\ney'J^{x,K^,X^,X2) 



9\ \2r^ i„,\/T>JJz 
27r2 

— ^^re//(a;, Ai, A2; x\ K^, A'l, A2) 



2(27r)^ 

x$4^^(x',/t'^,A;,A'2), (4.14) 



where Veff is the effective qq interactions in eq. (|2.19|) . 

Melosh transformation is exact only for free theory. With interactions incorporated, 
the use of Melosh transformation is only an approximation. This approximation may be 
reasonably good for the lowest spin bound states, such as the lowest-lying scalar and vector 
mesons, since the wavefunction (the valence quark amplitude in the light-front bound states) 
for them is a scale function to the rotational transformation. In other words, the "orbit" 
angular momentum which is dynamically dependent in the light-front formulation may not 
contribute to the total spin of these lowest spin hadrons. Moreover, in this paper, we focus 
on heavy quark systems. As we shall see next, due to the heavy quark spin symmetry, 
Melosh transformation results in the exact spin structure with fixed parity for scalar and 
vector heavy mesons. 
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B. Bound state equation for heavy quarkonia 



In this section, we shall explicitly consider the heavy quarkonium states. First of all, for 
quarkonia, the wave function ( |4.9| ) can be further simplified, especially for its spin structure 
due to the spin symmetry in HQET. Within the framework of light-front HQET, eq.( |4.9| ) in 
the heavy quark limit is reduced to: 



^<Poo(y^ ^^±M{h, Ai)4(^2, A2)). 



A1A2 



(4.15) 



Here the wavefunction (pQ-Q^y, may be mass dependent due to the kinetic energy in Hxo 
[see ( p.26|) ]. The Melosh transformation matrix element ( [4.11| ) in quarkonium states becomes 
a pure kinematic factor. 



pOO 

-"■A1A2 



2V2 



for a pseudoscalar meson, and 



(4.16) 



-"-AiAa 



u{v, Ai) ^{Jz)v{v, A2) 



for vector mesons. The light-front spinors for heavy quarks are given by 



so that 



u 



[v, X) = [1 + jwx 



Xa, 



v{v,X) = (l + 



a ■ v± — (3 



X-x, 



u{y, X)u{v, A') 



-^AA' 



Yu{v,X)u{v,X) 



1+ ^ 



v{v, X)v{v, A') = — ^5aa' , J2 ^(^' A)^^(^^, A) = ^, 
and the polarization vector is defined by 



(4.17) 



(4.18) 



(4.19) 
(4.20) 
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e^(±l) = (-|e^-i;^,0,ex) , e^(0) = -{d_l^y+ ^y^y ex(±l) = ± 0- (4.21) 

Thus we have constructed the hght-front heavy quarkonium bound states in the heavy mass 
hmit, which have definite spin and parity. The corresponding spin tensor structures are 
given by eqs. (|4.16| ) and (|4.17|) . 



Note that the heavy quarkonium states in heavy mass hmit are labelled by the residual 
center mass momentum K^^. We may normalize eq. ( |4.15| ) as follows: 



(vl>(K', J', J',)\^iK, J, J,)) = 2i27if K+5%K - K')5j,j5j,j^, (4.22) 
which leads to 



dyd K± 



\(l)g-g{y,K^)\' = l. (4.23) 



2(27r)3 

With the quarkonium states given above, it is easy to derive the corresponding bound 
state equation. In the weak-coupling scheme, Hlf = Hxo where Hxq is given by ( 3.251) . 
Thus, the quarkonium bound state equation in light-front HQET is given by 

{K--Hxo)\^{P,J,J,)) = 0. (4.24) 

The free energy part of the quarkonia states is extremely simple, 

K- - k- - k~ = -^2A' ^^(A'- Ml), (4.25) 

where Mq = —A is a residual invariant mass (=the invariant mass Mq subtracted by the 
mass dependent terms, here is just a notation rather than a real square of a quantity). 
It follows that eq.(|]2|) can be expressed explicitly by 



-2 



2A--A 



2 



2kI + k\2y' -2y + l) U^-^{y, k^) = [ - ^X'Cj Ine ) ^^^d/, k 



2(27r)3 [{y-y'y 

+ 7 ^^(^ - >^')]<t>QQ(y'^ -x)- (4-26) 

[K± - + {y- y'YK J 

This is the bound state equation for heavy quarkonia in the weak-coupling treatment of low 
energy QCD. 
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C. Bound state equation for heavy-light quark systems 



In the previous subsection, we have derived exphcitly the bound state equation for 
quarkonia. In the last few years, the heavy-light quark systems have been extensively ex- 
plored theoretically and experimentally. The discovery of the heavy quark symmetry in the 



heavy mass limit |rU| allows one to extract the Kobayashi-Maskawa matrix element \Vcb\ 
without knowing the detailed structure of the heavy-light mesons. However, to have a com- 
plete description for various heavy hadron processes, one has to know a number of heavy 
hadron matrix elements involved the heavy hadron bound states. Currently, most of heavy 
hadron matrix elements have only been calculated in various phenomenological models, such 
as quark models, QCD sum rules etc. It is necessary to find these heavy-light hadron bound 
states from the fundamental QCD. In this subsection, we shall derive from the low energy 
effective QCD Hamiltonian the bound state equation obeyed by the heavy mesons with one 
heavy quark. 

For the heavy mesons with one heavy quark, due to the heavy quark symmetry the general 
form of the wavefunction in the weak-coupling scheme can also be simplified. Considering 
the heavy quark mass limit, the wavefunction ( ^.9| ) can be expressed as 



xRiXiy^ ^±)<pQQiy^ ^Mik, Ai), d\p, A2)). (4.27) 



Here the Melosh transformation matrix element, eq.( [4.11| ), is a little more complicated in 
comparing to the heavy quarkonium state: 

1 



00 

A1A2 



for a pseudoscalar meson, and 



/, .iW(A), ^(^-^i)7^^(P>^^) (4-28) 



R 



Ai A2 



J^^l^^m, ^(J^>(P^ (4.29) 



2^2(yr + ^i±^) 
for vector mesons. The light-front spinors for the light antiquark is 
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' 1 ^ 



X-x, (4.30) 



and the polarization vector is still given by eq.( [4.21| ). 

Eq. (|4.27| ) is a light-front heavy-light meson bound state in the symmetry limit {rriQ — »• 
oo), which has the definite spin and parity. The corresponding bound state equation then 
becomes 

(A^ + (1 - ,)A^ - !i±^),.-.(,, i„ A„ A,) 



+ {K^fj ^l^vay-y\^^-^l)<^ii'{y',^l,>^u>^2), (4.3i) 

where Vqg is given by eq. ( |3.34| ), 



$^^^(7/, K^, Ai, A2) = 0Qg(y, K^)RiXiy, f^±), (4.32) 



and R\'l\2 is determined by eq. ( |4.28| ) or ( |4.29| ). Here the light antiquark is a brown muck, a 



current quark surround by infinite gluons and qq pairs that result a constituent quark mass 
which is a function of A. 

Thus, we have derived in this section the bound state equations in the weak-coupling 
scheme of the non-perturbative QCD for the light-light, heavy- heavy and heavy-light mesons. 
By solving these equations and comparing with experimental data, such as meson mass spec- 
troscopy, we can determine the A dependence of the constituent quark masses, the effective 
coupling constant as well as the wavefunction renormalization (anomalous dimensions) of 
hadronic states. Then we are able to use the corresponding wavefunctions to describe and 
predict various hadronic processes. The low energy cutoffs, or more precisely, the low energy 
scale dependences indeed reveal the inherent QCD dynamics of hadronic bound states. For 
completeness, we should also derive the bound state equation for glueball states {gg bound 
states). The glueball bound state equation is not only the basis for the study of the currently 
searching glueball states, but also allow us to determine another very important quantity in 
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the present weak-coupling treatment of low energy QCD, i.e., the constituent gluon mass and 
its scale dependence. But in this paper, we shall mainly consider heayy hadron systems. As 
we have seen, in the heavy quarkonium bound states, the constituent light quark and gluon 
masses do not appear. We only need to determine the effective coupling constant. Thus, the 
quarkonium states are the simplest systems in present theory. In fact, the determination of 
the scale dependence of the effective coupling constant g\ is the most important problem, 
from which we can test whether a weak-coupling treatment of nonperturbative QCD can 
be realized in this framework. In the following sections, we shall study heavy quarkonia in 
details. The extension to heavy-light quark systems will be briefly discussed and the more 
explicit study will be presented in the forthcoming publication. 

Before proceeding to solve the bound state equations derived in this section, we shall 
demonstrate first from these bound state equations how quark confinement is realized on 
the light-front. 

V. QUARK CONFINEMENT ON THE LIGHT-FRONT 

In the conventional picture, QCD has a complex vacuum that contains infinite quark pairs 
and gluons necessary for confinement and chiral symmetry breaking. On the light-front, the 
longitudinal momentum of physical particles is always positive, p+ — p'^+p^ > 0. As a result, 
only these constituents with zero longitudinal momentum (called zero modes) can occupy 
the light-front vacuum. The zero modes carry an extremely high light-front energy which 
has been integrated out in the similarity renormalization group scheme. Consequently, some 
equivalent effective interactions associated with the effect of the nontrivial QCD vacuum 
are generated in the low energy Hamiltonian Hx- Furthermore, the use of the constituent 
picture in the weak-coupling scheme forbids possible occurrence of any zero modes in Hx 
in the subsequent computations. Therefore, light-front QCD vacuum remains trivial. The 
nature of nontrivial QCD vacuum structure, the confinement as well as the chiral symmetry 
breaking must be made manifestly in Hx in terms of new effective interactions. In fact. 
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upon to the second order calculation in the similarity renormalization scheme (see sections 
II and III), the effective Hamiltonian Hx already contains a confining interaction. The 
interactions associated with the chiral symmetry breaking may be manifested in the fourth 



order computation of as pointed out by Wilson [|T5[, but these interactions are not 
important in the study of heavy hadrons. Hence, next we shall only discuss the quark 
confinement in terms of the light-front bound state equations, from which a light-front 
picture of confinement mechanism becomes transparent. 

To be specific, we take the following criteria as a definition of quark confinement: i) No 
color non-singlet bound states exist in nature, only color singlet states with finite masses 
can be produced and observed; ii) There is a confining potential for quark interaction such 
that quarks cannot be well-separated; iii) The conclusions of i-ii) are only true for QCD 
but not for QED. If conditions i-iii) could be verified from the low energy effective QCD 
Hamiltonian Hxq and the corresponding bound state equations, then quark confinement is 
realized on the light-front. Here we shall take heavy quarkonia as an explicit example. Some 
of the ideas for this confinement picture have been discussed in P,p^. 

In the present formulation of low-energy QCD, non-existence of color non-singlet bound 
states is essentially related to infrared divergences in the effective Hamiltonian. First of all, 
we shall show how only for physical states the infrared divergence in the quark self-energy 
correction is cancelled exactly by the same divergence from the uncancelled instantaneous 
interaction in eq.( |4.26|) . It is obvious that when y' — > the uncancelled instantaneous inter- 
action leads to a severe infrared divergence. Assuming that (pg-gly, k,±) is a smooth function 
with respect to y and k,±, and vanishes when y' ^ oo. Then the dominant contribution of 
the following integral is given by, 

= (|^(n(n In e) ct>QQ{yM). (5.1) 
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A more complete computation with an explicit light-front wavefunction will be given in the 
next section. 

Eq.( |5.1| ) indicates that in bound state equations, the uncancelled instantaneous interac- 
tion contains a logarithmic infrared divergence. Except for the color factor, this infrared 
divergence has the same form as the divergence in the self-energy correlation. From the 
bound state equation (^4.26|) , we immediately obtain the following conclusions. 

(a) . For a single (constituent) quark state, the bound state equation simply leads to 

]V^ = -|^C,ln. (5.2) 

This means that mass correction for single quark states is infinite (infrared divergent) and 
cannot be renormalized away in the spirit of gauge invariance. Equivalently speaking, single 
quark states carry an infinite mass and therefore they cannot be produced. 

(b) . For color non-singlet composite states, the color factor (T")q,/j(T")5^ in the QQ inter- 
action is different from the color factor dapCj = {T°'T"-)a(3- Therefore, the infrared divergence 
in the self-energy correction also cannot be cancelled by the corresponding divergence from 
the uncancelled instantaneous interaction. As a result, color non-singlet composite states 
are infinitely heavy so that they cannot be produced as well. 

(c) . For color singlet QQ states, the color factor (T")(T°) becomes {T°'T°') = Cf. Thus, 
the infrared divergences in eq. ([4.26|) are completely cancelled and the binding energy of the 



corresponding bound states is guaranteed to be finite. In other words, only color singlet 
composite particles are physically observable bound states, as a solution of eq. ([4.26|) . 

The above conclusion is also true for heavy-light and light-light hadronic states [see 
eqs.( [4.14| ) and ( |4.31|) ]. This provides the first condition for quark confinement in light- front 



QCD. Indeed, the physical origin of the above result is very clear. Light-front infrared 
divergences are associated with violation of gauge invariance. Only in gauge noninvariant 
sectors, light-front infrared divergences may occur. In gauge invariant sectors, infrared 
divergences must be automatically cancelled. Therefore, the above conclusion is a natural 
consequence of gauge invariance. 
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Physically, in order to be consistent with the above conclusion, confinement must also 
imply the existence of a confining potential so that quarks cannot be well-separated to 
become asymptotically free states. Now we can show that the interactions in effective 
Hamiltonian ^3. 2^ contains indeed both a confining and a Coulomb potentials. 



The Coulomb potential can be easily obtained by applying the Fourier transformation to 
the second term in (|3.27|) . It is convenient to perform the calculation in the frame K± = 0, 
in which 



i / <;.-<iVe"-— ' (A) ^ (5.3) 



x=! (x- f 

where = — k'^ = K^{y — y'), q± = ki±_ — /c3j_ = n± — for K± = 0. Eq.( ^73|) shows 
that the Coulomb potential on the light-front for quarkonium states has the form 

= -fcf^ ^ = -|ic^Ai, (5.4) 

where 



n = yxi + (A)^(^-)2 (5.5) 

is defined to be a "radial" variable in light-front space [0. 

The confining potential corresponds to the finite part of the non-cancelled instantaneous 
interaction in ( p.27| ). Its Fourier transformation is relatively complicated. The general 
expression is 



27r2 7o \xi_\q^ra 



where q±_rn = \/ ■^q'^ ~ Y+^q'^'^, and Ji{x) is a Bessel function. An analytic solution to the 
integral (|5.6|) may be difficult to carry out. However, the nature of confining interactions is 
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a large distance QCD dynamics. We may consider the integral for large x~ and x±. In this 
case, if q^x~ and/or are large, the integration vanishes, yet Ji{x) = | + + ■ ■ ■ for 

small X. The dominant contribution of the integral ( p.6| ) for large x~ and x± comes from 
the small such that q^x" and/or must remain small and therefore 



This corresponds to 



^.,+..-2Ji(k±k±J (5.7) 
\x±\q±m 



, 1 , , K+ , , 

q < — and/or q < — -r. (5.8) 

x~ \x±rX^ 



If g+ < ^ < i^jJx?, eq.dU) is reduced to 



'Kc, <i,+-L(il,+ - A_,«) =-4^C,(l„|x-| + ln.). (5^9) 



-2 

where a term ~ ^ is neglected since x~ is large, and e is an infrared cutoff on g+. The 
infrared logarithmic divergence (~ Ine) exactly cancels the divergence in the self-energy 
corrections in H\, so that the remaining is a logarithmic confining potential: 

Kon/.(a;-,x^) = ^^Cf\n\x-\. (5.10) 

Similarly, when g+ < ^^^^2^2 < have 

where the term ~ ^ has also been ignored because of large x\. Again, the infrared diver- 
gence (~ Ine) is cancelled in Hx, and we obtain the following confining potential: 

ycon,{x-,x,) = ^C,\n^. (5.12) 

Thus, the effective Hamiltonian Hxo contains a logarithmic confining potential in all the 
directions of x~ and Xj_ coordinates. Note that for heavy quarkonia, a logarithmic confining 
potential provides indeed a good description to the spectroscopy and leptonic decays . 



More details of computation will be given in the next section. Nevertheless, we have explic- 
itly shown here that if^o exhibits a Coulomb potential at short distances and a confining 
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potential at long distances. The second condition for quark confinement is verified on the 
light-front. 

Finally, we shall argue that the above mechanism of quark confinement is indeed only 
true for QCD. As we have seen the light-front confinement potential is just an effect of the 
non-cancellation between instantaneous interaction and one transverse gluon interaction. 
Such a non-cancellation arises from the requirement in the similarity renormalization group 
scheme that the transverse gluon energy cannot be below a certain value (about a few 
MeV). This requirement is naturally satisfied if the gluon mass is nonzero in low energy 
scale. Unlike the constituent quark mass which we know is an effect of the spontaneous 
chiral symmetry breaking, the origin of constituent massive gluons is not very clear at 
present. The assumption of massive gluons here may also violate gauge invariance but it is 
not unnatural. In fact, if gluon were massless like photons, the hadronic spectra would be 
continuous rather than discrete, as Wilson pointed out recently [|15|. A typical evidence of 
gluons being massive in the low energy domain is the possible existence of glueball states 
which is still a very active topic in current experimental searches P3[. The massive gluon 



must be originated from the nonlinear interactions in non-abelian gauge theory. Therefore, 
the non-cancellation of the instantaneous interaction in the low energy domain is indeed 
a consequence of the existence of the constituent massive gluons due to the non-abelian 
gauge interactions. This is independent of any particular renormalization scheme. The use 
of the low energy cutoff A just gives us a simple realization of this confining picture that 
the massive gluon exchange energy cannot run down to the zero value in nonperturbative 
QCD. In the case of lacking the mechanism of how the massive gluons are generated, the 
determination of the gluon mass lies on the solution of bound state equations. 

Based on the above discussion, it is now easy to find that the confinement mechanism 
described in this section is not valid in QED. First of all, the infrared divergence in the 
self energy is also a result of the noncancellation between the instantaneous interaction 
self-energy diagram and the one-loop self-energy corrections [see eq.( ^.28D ]. In QED, since 



the photon mass is always zero, the photon energy in the one-loop self-energy correction 
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covers the entire range from zero to infinity. Thus, in QED, we can always choose the low 
energy cutoff A being zero. (We shall further explain in Section VII that this is indeed the 
only choice for applying similarity renormalization group approach to QED. Otherwise the 
resulting effective QED theory is inconsistent with the perturbative QED theory.) Then a 
direct calculation for ( |3.28| ) with A = shows that the infrared divergences do not occur 
in the electron self-energy correction. As a result, the renormalized single electron mass 
is finite, in contrast to the divergent mass of single quark states. For the same reason, 
with A = 0, the instantaneous interaction in the effective QED Hamiltonian is also exactly 
cancelled by the same interaction from one transverse photon exchange so that only one 
photon exchange Coulomb interaction remains. Therefore, using similarity renormalization 
group approach to QED, we obtain a conventional effective QED Hamiltonian which only 
contains the Coulomb interaction. Such an effective Hamiltonian is the basis in the study 
of positronium bound states. More discussion will be given in Section VII. 

Now we shall study how a weak-coupling treatment scheme works in solving hadronic 
bound state problem in the present low energy QCD formulation. 



VI. QCD DESCRIPTION OF QUARKONIA ON THE LIGHT-FRONT 



A numerical computation to the heavy hadron bound state equations, eqs.( [4.26"D and 



( [1.31| ), is actually not too difficult. However, to have a deeper insight about the internal 
structure of light-front bound states and to determine the scale dependence of the effective 
coupling constant in H\, it is better to have an analytic analysis. In this paper, the light- front 
wavefunction ansatz will be used to solve the bound state equations for heavy quarkonia, 
from which some general properties of the low energy scaling in the similarity renormalization 
group can be extracted. It also provides a direct test whether the weak-coupling treatment 
of nonperturbative QCD can be realized. 
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A. A general analysis of light-front wavefunctions 



For heavy quark systems, the wavefunctions considered in the previous section are de- 
fined in heavy mass hmit. Most of the l/mg corrections can be handled in the standard 
perturbation theory in the present framework, except for the kinetic energy for quarkonia. 
The heavy hadronic wavefunctions in the heavy mass hmit can be tremendously simplified. 

First of all, the heavy quark kinematics have already added some constraints on the 
general form of the light-front wavefunction k,±). The kinetic energy part (the left hand 
side of these bound state equations in section IV) shows that when we introduce the residual 
longitudinal momentum fraction y for heavy quarks, the longitudinal momentum fraction 
dependence in is quite different for the heavy-heavy, heavy-light and light-light mesons. 

For the light-light mesons, such as pions, rhos, kaons etc., the wavefunction (t)qq{x,n±) 
must vanish at the endpoint a; = or 1. This can be seen from the kinetic energy contribution 
in the bound state equation [see eq. ( [4.14[ )], 



X I — X 

To ensure that the bound state equation is well defined in the entire range of momentum 
space, \(f)qq{x, K±)\'^ must fall down to zero in the longitudinal direction not slower than 
1/x and 1/(1 — x) when x — > and 1, respectively. In other words, at least 0gg(a;, k^) ~ 



x{l — X) . 

For heavy-light quark mesons, namely the B and D mesons, the wavefunction (pQqiy, k,±) 
is required to vanish aty = 0, where y is the residual longitudinal momentum fraction carried 
by the light quark. This is because the kinetic energy in eq. ( |4.31| ) contains a singularity at 

y = o, 

- Ml - = A% (1 - y)-K^ - ^i±Z!i. (6.2) 

y 

On the other hand, since < ?/ < cxo, (pQ-qiy, k_l) should also vanish when y oo. Hence, a 
possible simple form is (pQq{y,K±) ~ i/ye""^ or y/ye~"^^ The y dependence in 4>Qq{y,K±) 
is obviously different from the x dependence in (pqg{x, k±). 
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For heavy quarkonia, the kinetic energy in the corresponding bound state equation (f4.26|) 

is: 



M2 - Ml ^K^-mI = 2K'- — 



2k\ + A (2y2 -2y + l) . (6.3) 

rriQ L J 

Since — oo < y < oo, the normahzation forces (pQ-qiy^ to vanish as y ^ ±00. Therefore a 
possible form is (pQQ^y, k±) ~ e~"^^ Obviously, in heavy quark mass limit, the y dependence 
in 4>QQ{y, K±) is very different from the above two cases. 

On the other hand, the transverse momentum dependence in these light-front wavefunc- 
tions should be more or less similar. They all vanish at k,± ±00. A simple form of the 
K_L dependence for these wavefunctions is a Gaussian function: e~'^J-/^'^^. 

The above analysis of light-front wavefunctions is only based on the kinetic energy prop- 
erties of the constituents. Currently, many investigations on the hadronic structures use 
phenomenological light-front wavefunctions. One of such phenomenological wavefunctions 
that has been widely used in the study of heavy hadron structure is the so-called BSW 
wavefunction, introduced by Bauer et al. p4| . 



^Bsw{x, = M\/x{l - x) exp ( I exp 



2(3; - ^oY 



(6.4) 



2a; 

where Af is a normalization constant, u a parameter of order Aqcd, = (| — "m^^ ) ' 
and Mh, mi, and m2 are the hadron, quark, and antiquark masses respectively. In the 
phenomenological description, the parameters uj, and {i = 1, 2) in (|6.4|) are fitted from 
data. Here we are of course interested in the dynamical determination of these parameters. 

As we have pointed out in passing, for heavy quark systems, the l/mq corrections can 
be well treated perturbatively in our framework (except for the kinetic energy of quarkonia). 
Here we are only interested in the solution of the wavefunctions in the heavy mass limit, 
where eq. ( |6.4| ) can be further simplified. 

Explicitly, for heavy-light quark systems, such as the B and D mesons, one can easily 
find that in the heavy mass limit, 

mi = TTiQ ~ Mh, iTiq « mq so that Xq = 0. (6.5) 
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Meanwhile, from eg. ( |3.31| ), we also have 

Mhx = Ay. (6.6) 



Furthermore, the factor J x{l — x) can be rewritten by y/y in according to the discussion 



on eq. (|6y^ ). Thus, the BSW wavefunction is reduced to 



(f>Qqiy, t^±) = ^fy/y exp y-^j j ■ ^^■^'> 

This agrees with our qualitative analysis given before. Using such a wavefunction we have 
already computed the universal Isgur-Wise function in _B ^ D,D* decays 

i{v ■ v') = (6.8) 
V ■ v' 

and from which we obtained the slope of ^(f ■ v') at the zero-recoil point, = — ■^'(1) = 1, 
in excellent agreement with the recent CLCO result [^] of = 1.01 ± 0.15 ± 0.09. This 
result is independent of the value of u, and therefore is independent of further dynamics of 
QCD involved in the corresponding bound state equation. The simple form ( |6.7|) is just a 
consequence of heavy quark symmetry in the our low energy theory. We may argue that 
= 1 could be a universal identity. 

For heavy quarkonia, such as the bb and cc states, mi = 777,2 = itT'Q which leads to = 1/2 
in eq.( |6.4| ). Also, the longitudinal momentum fraction in ( |6.4| ) is defined by x = Pi /P^, its 
relation to the residual longitudinal momentum fraction is given by 

MH{x-^) = Ay. (6.9) 



In addition, the factor yx{l — x) must be totally dropped as we have seen from the discussion 
on eq. (|6.3| ). Therefore the BSW wavefunction for quarkonia is reduced to 



,2 



A^ 



_ iy, «:x) = AT exp [-^j exp \^-^y'j , (6-10) 

which is the exact form as we expected from the qualitative analysis. Here we have not taken 
the limit of tuq 00 for heavy quarkonia. Thus a possible mq dependence in wavefunction 
may be hidden in the parameter u. 
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Another phenomenological light-front wavefunction which has been widely used for both 
light and heavy mesons has the form [EH] 



V,„(.,.J.W|iexp(-|l) exp(-^], (6.11) 



where is defined by 



x=^i±^, l-x='-^-^; e, = ^Kl + Kl + m^ (z = l,2) (6.12) 
ei + 62 61 + 62 ^ 



as a pretended ^-component of relative momentum while ydKz/dz is the Jacobian of trans- 
formation from (x, to k = {k,±,Kz)- This wavefunction has been used frequently in 
various studies of hadronic transitions. In particular, it has been shown that this wavefunc- 
tion describes satisfactorily the pion elastic form factor up to ~ 10 GeV^ . 
For heavy quarkonia with mi = m2 = mg, we may have 

= Mo(x-i) ^Ay. (6.13) 



Thus, ydKz/dx = constant, and eq.( |6.11| ) is reduced to the same form obtained from the 
BSW wavefunction, 

(j)gg{x, K±) — > (f^Q-qiy, = A/" exp exp ^-^Z/^j • (6-14) 

Therefore, the above Gaussian-type ansatz should be a very good candidate for the low-lying 
quarkonium states. In the following, we start with this wavefunction ansatz to solve the light- 
front quarkonium bound state equation, and from which to determine the low energy scaling 
dynamics and develop the weak-coupling treatment of the heavy hadron bound states. 

B. A weak-coupling realization of the nonperturbative QCD description for heavy 

quarkonia 

Based on the analysis in the last section, we take the normalized wavefunction ansatz of 
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0QQ(2/,/t±) = 4VA 



7r 



3/4 



exp 



K 



exp 



(6.15) 



as a solution (a trial wavefunction) of the heavy quarkonium bound state equation ( [4.26[ ). 
Note that here we have also specified the scale dependence of the wavefunction through 
the scale dependence of the parameter u!\. Substituting the above wavefunction into the 
quarkonium bound state equation ( [4 .261) and introducing the new variables 



Z = ^{y + y'), z = y-y', 



(6.16) 



we have 



^ I 2(27r)3 



exp 



X J e{X'\z\-qi-z'A^) + 



e{qi + z'k-XM) 



kin 



^^A^C/ In e + Enonc + ^CouU 



(6.17) 



where Skin represents the kinetic energy, Snonc is the contribution of the noncancellation of 
the instantaneous interaction, and Scoui from the Coulomb interaction. Furthermore, it is 
not very difficult to compute that 



dz(Pqj_ 1 



1 



2x2^ 



2tt' 



2(27r)3 z2 

:CfX^L + In ^ + Ei{zu^) + ^Erf(ti7)|, 



{qi + z'K 



(6.18) 



Coul 



dzcPqx^ A 



-^9lCf / ^--^^(gi + z'T - XM) exp {-^M + ^'A^) 



2(27r)3 ^2 + ^2A^ 



9l 



27r2 







1 r 


[VI 


1- Erf(ro)' 


+ — 


\ w 







(6.19) 



where 7 = 0.57721566... is the Euler constant, e is the small longitudinal momentum cutoff, 
the dimensionless w is defined by 
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(6.20) 



and El and Erf are the exponential integral function and the error function, respectively, 

2 i-x 



Ei(x) = / — dt , Erf(x) = ^ / e"* 
Jx t v^r JO 



(6.21) 



We may rewrite the term In j% in Enonc as 



, ^'e , , 2 , a' 
m -— J = In e + m +ln— . 
4a;| A2 



(6.22) 



It shows that Snonc contains a logarithmic divergence Ine which exactly cancels the same 
divergence from the self-energy correction, as expected, and the term In ro^ is the logarithmic 
confining energy. 

After the cancellation of the infrared Ine divergences in eq.( |6.17| ), the binding energy for 
heavy quarkonia is given by the kinetic energy plus the potential energy: 



2A =£, 



kin ~t~ ^conf ~l~ ^Coul 

A 



— {'iujl + h') + ^CfX'l 7 + Intu' + In ^ + E^{w') + ^Erf(tu 



9l 



a' 



27r' 



zu 



27r2 



9l 





1 - Erf (u;)' 


1 

+ — 


1 - e-'' 











a' 



2^-C,A^ F(t^) + ln ^ 



(6.23) 



where 



zu 



1 - 2Erf (t<7) 



1 



1 -e" 



(6.24) 



which is a dimensionless function. In Fig.l, we plot the confining potential energy, the 
Coulomb potential energy and the totally potential energy as functions of zu which is pro- 
portional to the radial variable in light-front space. 



zu ~ — ~ ri. 



(6.25) 



Fig.l shows that the total potential energy is a typical combination of the Coulomb potential 
in short distance and a confining potential in long distance that has been widely used in 
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previous phenomenological describing hadronic states, but it is now explicitly derived from 
QCD. Furthermore, eq.( |6.23| ) also indicates that without considering the kinetic energy, we 
cannot find stable quarkonium bound states. The kinetic energy balances the potential 
energy and ensures the existence of a stable solution for ( |6.23| ). Therefore, the first order 
kinetic energy in HQET is an important nonperturbative effect in binding two heavy quarks. 



as noticed first by Mannel et al. in their attempt of applying HQET to heavy quarkonium 
system. 

If we know the experimental value of the quarkonium binding energy A, minimizing 



eq. (|6.23|) can completely determine the parameter ux and the coupling constant gx. The 
precise value of quarkonium binding energy that can be compared with the data in Particle 
Data Group must include the spin-splitting energy {l/mq corrections) which we will 



present in the forthcoming paper |2^. Here, to justify whether a weak-coupling treatment 



of the nonperturbative QCD can become possible in the present formulation, we will give 
a schematic calculation. It is known that A is of the same order as Aqcd which is about 



100 ~ 400 MeV. To solve (|6.23|) we shall take several values of A within the above range. 



We choose the low energy cutoff about a typical hadronic energy, A = 1 GeV. The charmed 
and bottom quark masses used here are = 1.4 GeV and = 4.8 GeV. The results are 
listed in Tables I and II for charmonium and bottomonium, respectively, where uxo denotes 
the minimum point of the binding energy ( |6.23|) . 

We see from the Tables I-II that the coupling constant ax = g'x/'^T^ is very small. For 
instance, with A = 200 MeV, we obtain 

{ 0.02665 charmonium, 
(6.26) 
0.06795 bottomonium, 

which is much smaller than that extrapolated from the running coupling constant in the 
naive perturbative QCD calculation. The parameter uoxo in the quarkonium wavefunction 
is the mean value of the (transverse) momentum square of heavy quark inside the heavy 
quarkonia: 
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{kl) = ool (6.27) 

For charmonium, we can see that the resulting ojx are typical values of Aqcd ~ A. The 
kinetic energy is about a half of the potential energy. For bottomonium, we find that the 
binding energy A cannot be too large. In fact, when A is over about 260 MeV, eq.( |6.23D 



has no solution. Meanwhile, compared to charmonium, the effective coupling constant is 
relatively large (in contrast to the perturbative running coupling constant which is smaller 
with increasing mg if it is taken as the mass scale). Also the values of ux in bottomonium 
wavefunctions are larger than that in charmonium. The difference between charmonium 
and bottomonium in the nonperturbative calculation may be understood as follows. As we 
know, in the nonrelativistic quark model, the quark momentum in quarkonia is proportional 
to the quark mass, uj\ ~ mq [^. Our relativistic QCD bound state solution exhibits such 
a property. This is why the values of ux for bottomonium are much larger than that for 
charmonium. As a result, the bottomonium kinetic energy (~ uD becomes large as well. 
To have a nonperturbative balance between the kinetic energy and the potential energy 
in the bound states, the coupling constant in bottomonium must be larger than that in 
charmonium. All these properties now have been manifested in the solution of eq. ( |6.23| ). A 
more precise determination of ax (i.e., gx) requires an accurate computation of the low-lying 
quarkonium spectroscopy with the l/mg corrections included PS|. Nevertheless, it has been 



shown that the effective coupling constant in the low energy Hamiltonian Hx is very small 
at the hadronic mass scale. 

In order to see how this weak coupling constant varies with the cutoff A, we take A = 200 
MeV and vary the value of A around 1 GeV. The result is listed in Table III. We find that 
the coupling constant is decreased very faster with increasing A. In other words, with a 
suitable choice of the low energy cutoff A in the similarity renormalization group scheme, 
we can make the effective coupling constant ax in Hx arbitrarily small, and therefore the 
weak-coupling treatment of the non-perturbative QCD can be achieved in terms of Hx such 
that the corrections from Hxi can be computed perturbatively. Thus, we have provided the 
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first explicit realization of recently proposed the weak-coupling treatment of nonperturbative 
QCD on the light-front [g. 

We must emphasize here that a\ is not the physical coupling constant = (7^/471. The 
later is of order unity at the hadronic mass scale. A detailed analysis of the A-dependence 
and the relation between ax and will be discussed in the next. 

VII. SIMILARITY RENORMALIZATION GROUP EQUATION AND LOW 
ENERGY RUNNING COUPLING CONSTANT 

In this section, we shall discuss the scale dependence of the coupling constant, the con- 
stituent quark and gluon masses as well as the wavefunctions. For heavy quarkonia, the 
bound state equation does not include the constituent gluon mass. The heavy quark mass is 
larger than the usual hadronic mass scale. Its A-dependence should be very weak that can 
be neglected in the present discussion. Thus the remainings are the A-dependence of the 
coupling constant gx and the wavefunctions, the later is described through the A-dependence 
of the parameter ux- 

From these solutions in Tables I to III, we find that the values of dimensionless parameter 
zu = are greater than 2.5. When x > 2.5, the exponential integral function and the 

error function are simply reduced to Ei(a;) = and Erf(a;) = 1. Thus, the dimensionless 
function F{w) can be expressed approximately by 




zu > 2.5, 



(7.1) 



with an error less than 10 ^. Hence we can simply rewrite eq.( |6.23| ) as 




(7.2) 



Minimizing A with respect to ux, we obtain 




(7.3) 



Therefore, eq. (|7.2| ) becomes 
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where uxo is a solution of ( [7.3|) . Eqs.(|7l3|) and ( |7.4|) determine the A-dependences of the 
couphng constant gx and the wavefunction parameter ujxq- 

Directly and analytically solving eqs.ffOI) and (|7^) is not obviously possible. The 
nonperturbative balance between the kinetic energy and the potential energy implies that 



LJx ~ ynigA. Meanwhile, since it is the binding energy of heavy quarkonia, A should be 
A-independent. We then obtain 





TT 






TT 









'A^ / A 



AVV 2mQ;i+^ + in^-lnc^2^ + 0F^ 



1 



a + 61n^ 



(7.5) 

A"' 

where the coefficients a and b can be obtained by numerically solving eqs. (|7.3|) and ( [7.4|) . 
The coefficient b is almost a constant (with a slight dependence on mq but independence 
on A and A), while a depends on A, A and also rriQ. For A > 0.6 GeV, the A-dependence of 
the parameter a is negligible. In Fig. 2, we plot the A-dependence of the effective coupling 
constant ax for charmonium. The dots are the numerical solutions of ( |6.23|) and the solid 



line is given by the analytical form (7^) with b = 1.15, and a = —0.25 for A = 0.2 GeV and 
a = 1.1 for A = 0.4 GeV. We can see that ( |7.5| ) is a very good analytical solution of the 
eqs.([7.3[) and ( [7.4[ ) [or of the minimizing eq.( |6.23| )]. 

The above solution shows that with increasing A, ax becomes weaker and weaker. Mean- 
while, we also find that the confining energy becomes more and more dominant in the binding 
energy (See Table IV). Fig.3 is a plot of the parameter ujxo as a function of A, from which 
we also see that with increasing A, uxo is decreased. Correspondingly, the distance between 
two quarks inside the quarkonia, r/ ~ is increased. This is why the confining interaction 
becomes more and more important. On the other hand, the confining interaction comes from 
the noncancellation of instantaneous gluon exchange with energy below the scale A. With 
the larger A, the more the instantaneous interaction contributes to Hx- Thus, the above 
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conclusion can also be directly understood from the low energy Hamiltonian H\. Compared 
to the canonical QCD theory, the confining interaction should become more important if 
the scale would be smaller, and correspondingly the running coupling constant 
becomes larger. This indicates that there is an inverse correspondence between the effective 
coupling constant ax in the low energy Hamiltonian Hx and the running coupling constant 
a(Q2) in the full QCD theory: 



In other words, the weak- coupling treatment of the low energy confining Hamiltonian Hx may 
correspond to an inverse strong- coupling expansion of the full QCD theory. The similarity 
renormalization group approach provides an implicative realization for such an expansion. 
This may be the inherent property why the nonperturbative QCD can be treated as weak- 
coupling problem in the similarity renormalization group scheme and why we can find the 
confining interaction in this weak-coupling QCD formulation. 

We also find from Table IV that the confining interaction plays a more important role 
than the Coulomb interaction in the determination of the quarkonium bound states. This re- 
sult is different from the usual understanding in the nonrelativistic phenomenological descrip- 
tion that the dominant contribution in heavy quarkonium spectroscopy is the Coulomb inter- 
action. This discrepancy can be understood as follows. The currently relativistic light-front 
description for heavy quark system mostly uses the heavy quark masses of mc = 1.3 ~ 1.4 
GeV and nib = 4.8 GeV or less (In Particle Data Group [^, rrii, = 1.0 ~ 1.6 GeV and 
rric = 4.1 ~ 4.5 GeV). Thus, the heavy quarkonium binding energies, A = Mh — 2mQ, 
might be positive [the lowest charmonium ground state M{ric{lS)) = 2.98 GeV, and the 
bottomonium M(T(1S')) = 9.46 GeV]. Therefore, the Coulomb energy is obviously not im- 
portant. The dominant contribution for binding quarkonium states must come from the 
nonperturbative balance between the kinetic energy and the confining energy. While, in the 
nonrelativistic phenomenological description, one used the larger quark masses, rric > 1.8 
GeV and rrih > 5.1 GeV |2^, such that the binding energy is negative and therefore the 



1 



and ~ — . 



(7.6) 
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Coulomb interactions must be dominant in this picture. Of course, on the hght-front, the 
structure of the bound state equation is different from the nonrelativistic Schrodinger equa- 
tion. There is no direct comparison. A real solution to the above discrepancy may be 
obtained after including the spin-splitting interactions (1/mQ-corrections). 

Now we can study the running behavior of the coupling constant in the similarity renor- 
malization group scheme. Denote 

A = A{gx,uJx,X). (7.7) 

The invariance of the binding energy A under the similarity renormalization group transfor- 
mation means that the A determined from Hx and Hy must be the same for A 7^ A'. Let 
A' = A + 6X, we obtain the corresponding similarity renormalization group equation 

where the quantity (3 is the similarity renormalization group j3 function which is defined by 



(7.9) 

A=A(gA) 



and is an anomalous dimension that describes the running properties of the bound state 
wavefunction. The /3 function can be computed from eq. ( |7.5| ), 

2b 



P = -9x I 



a + 2&ln^ 



A=A(gA) 

^ —gx (for a relatively large A >> A). (7-10) 
On the other hand, the running coupling constant in full QCD theory is given by 

t=r^ (7 11) 

where t = Mn%, and is a space-like momentum (the same as A^). Since the similarity 



renormalization group jS function of eq.( |7.10| ) is determined in the physical sector of low 
energy QCD dynamics, the low energy (3 function of the running coupling constant g{Q^) 
in the full theory should behave qualitatively the same. With this assumption, the /3{g) 
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function in the above equation may take the same form as eq.( [7.10|) in low momentum 
transfer (namely « Ag^^^^). This leads to 

fm=9l^, , Q'«^%CD. (7.12) 

and gl = gli^n'^) is a fixed coupling constant at the hadronic mass scale [j?. In terms of the 
running coupling constant a = g'^/Air, we have 

«(Q^) = «.(/x^)^ = co^, (7.13) 

where cq = as{^i^)^i^ / Is^^cd- This is consistent with eq . (17.61) . Furthermore, we see that 
the fixed point of the coupling constant under renormalization group transformation is the 
origin of g, and it is an infrared unstable (UV stable) fixed point, in consistence with the 
asymptotic freedom of QCD. 

To give a qualitative determination of the coefficient cq, we consider A = 0.75 ~ 1.5 GeV 
and A = 0.2 GeV, then = (14 ~ 56)A >> A . Rewriting ( |7.5| ) as the same form of 
( [7.13| ), we obtain: 

-r-2 

«A = (1.0-1.5) —. (7.14) 

-r2 



The corresponding ~ << A ~ ^^cd- From (|7.6|) , We may require that 



"A Q 



2 



(7.15) 



a(g2) \2 ■ 

It follows that Co = 1.0 ~ 1.5, namely for Q"^ « Ag^^ 

a(Q2) = (1.0~1.5)^. (7.16) 

This is just a qualitative estimation of the running coupling constant in the full QCD theory 
in low momentum transfer. A more accurate result may be obtained by exactly solving the 
/^-function of eq.( |7.10| ). The running coupling constant in high momentum transfer is given 
in the usual perturbative QCD theory. A light-front perturbative QCD calculation of the 
leading order running coupling constant can be found from Ref . [jl^ . The result is standard: 
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2N as(/U^) 12vr 



"^^^ l + 6o«s(/i2) In gV/i2 (33-2iV^)lngVA|^^- ^^'^^^ 

Up to date, no one precisely knows how the QCD couphng constant varies in low energy 
scale. However, it is interesting to see that the running coupling constant given by eqs.( [r.l3"D 



and ( [7.171 ) for the small and large respectively is indeed the basic assumption of the 
Richardson QQ potential |]30| : 

Vm = -Cf^. (7.18) 

where 

"^^'^ = (33-2iV,)ln(l + gVA^^z.)- ^^'^^^ 
The Richardson QQ potential is proposed to exhibit the asymptotic freedom of QCD in 
short distance and a linear potential in large distance. From eq. (|7.18| ), we see that for large 
g2 (g2 » A^^^), 



127r 

(33-2iv;yhrg7A|^^' 



which reproduces the result of the asymptotic freedom of QCD. The corresponding potential 
is just the Coulomb potential. For small Q^ (g^ << Aq,^^), 



and the corresponding potential from ( [7.18 ) becomes 



_ l27,C,klcD 1 



which is a Fourier transformation of the linear potential, 

V{r) = kr. (7.23) 
The Fourier transformation of ( [7.18| ) is the QQ potential in coordinate space: 
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roo dq 



-It 



where f{t) = [l-AJ^ yWUF^W^ 
used to describe quarkonium dynamics. 



pO| . The Richardson potential has successfully been 



Comparing with eqs.( |7.13D and ( [7.21|) , we have from the Richardson QQ potential (with 

12% 



Co 



1.4 



(7.25) 



(33 - 2Nf) 

This result agrees very well with eq.( |7.16| ). Consequently, although it is a very rough qual- 
itative analysis, the above result may imply that the confining Hamiltonian derived from 
light-front HQET in the similarity renormalization group scheme has also covered the dy- 
namics of linear confining potential. As we have known, phenomenological potential quark 
models based on the Richardson potential, the linear plus Coulomb potential (also called the 
Cornell potential) as well as the logarithmic potential all give a good description of quarko- 
nium dynamics Hence, it should be not surprising if our QCD confining Hamiltonian 
encompass the dynamic behavior of all these potentials. 



VIII. DISCUSSIONS ON THE FULL QCD VIA EFFECTIVE THEORY 

Thus far, the main ideas of the weak-coupling treatment on nonperturbative QCD pro- 
posed in the recent publication have, at least qualitatively, been achieved for heavy 
quarkonium. The low energy nonperturbative QCD theory is defined by the effective low 
energy Hamiltonian Hx. The key to solve this theory is to determine from the bound state 
equation the A-dependence of the effective coupling constant ax. Our result indicates that 
the low energy effective QCD Hamiltonian exhibits an alternative realization of the inverse 
strong coupling expansion of the full QCD theory. Thus, with a suitable choice of the cutoff 
A, the effective coupling constant ax (as well as gx) can be arbitrarily small. As an example, 
one may take A to be a constituent gluon mass (about a half of the glueball masses, such as 
the recent possible evidences of /o(1500) and ^(2230) P5|), 

A : (0.75 ~ 1.5) GeV. (8.1) 
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Then the effective couphng constant (with A = 200 MeV) is 

ae(A2) = 0.06 ~ 0.01 (8.2) 

which is very smalL The residual interaction Hxi is expanded in terms of this small coupling 
constant so that the corrections from H\i can be perturbatively computed, and the weak- 
coupling treatment of nonperturbative QCD is explicitly realized. 

Now the question is in what limit the effective theory can return back to the full theory 
of QCD so we can ensure that the present formulation is a complete consistent theory 
describing low energy QCD dynamics. This question is also important in the sense that 
with the cutoffs being introduced and the assumption of gluon being massive, the gauge 
symmetry and rotational symmetry may be broken down in the effective theory. Then we 
must know when all these symmetries can be restored. In ref. we have argued that when 
Qx Qsi the effective theory must recover the full QCD dynamics so that all symmetries 
are restored as well. However, it is not clear how this limit can be achieved. Here we shall 
attempt to answer this question from the heavy quarkonium solution. 

We have used three cutoffs in this paper: the UV cutoff A, the low energy cutoff A and the 
infrared longitudinal momentum cutoff e. The UV cutoff is renormalized away in ordinary 
perturbation theory so that there is no any explicit A dependence in our formulation. The 
longitudinal infrared cutoff is automatically cancelled in all the physical sectors, due mainly 
to the gauge invariance as we have seen from the calculations throughout this work. The 
final formulation only contains the low energy cutoff A. This A dependence is essentially 
associated with nonperturbative QCD dynamics. However, the similarity renormalization 
group invariance on physical observations allows us to further remove away the A dependence 
in physical sectors. We may define that the naive A — > limit brings the effective theory 
back to the full theory of QCD. 

How can the limit of A ^ theoretically bring the effective theory back to the full QCD 
theory? Firstly, recall that introducing the low energy cutoff in the effective theory is based 
on the assumption that gluon is massive in the low energy domain. Then the limit of A — is 
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only allowed when the gluon mass goes to zero. With this limit the broken gauge symmetry 
due to the massive gluons is now restored. Secondly, once the cutoff A is renormalized 
away, there is no any explicit cutoff dependence in the theory. Therefore the broken Lorentz 
symmetry due to the use of the explicit cutoff must be restored as well. Furthermore, the 
limit X —* corresponds to — > oo, where all gluons and quarks become current ones. 
Once all symmetries are restored and the current picture reemerges, the resulting theory 
should be the full QCD theory. 

We may first check what happens if we perform the same procedure to QED for positro- 
nium. With A — > 0, the similarity renormalization approach leads to an effective QED 
Hamiltonian in which the nonperturbative part Hxo only contains the Coulomb interaction, 
and the remaining is the radiative correction Hxi- No confining interactions and no infrared 
divergences occur, as expected. This is just the full QED theory used in the description of 
QED bound states. We can explicitly examine the above conclusion from eg. ( |6.23|) . To do 
so, we may first assume that the QED coupling constant is almost A-independent because 
it is always very weak in the whole range of energy scale. Thus, with A — 0, we have from 
eq.(S), 

_ 

£nonc. = 0, Scoul. = :TCf^/7^UJA, (8.3) 

TT 

namely, only Coulomb force contributes, and the confining interaction disappears, while the 
infrared divergence in the self-energy correction does not occur for A = 0. Combining with 
the kinetic energy (where the term A /uiq in £kin is very small so that it can be neglected), 
the totally binding energy is given by 



A = -^^2 - f-.CfV^iv. (8.4) 



Minimizing A with respect of u, we obtain 

^o--^Cf, A--^^Cf. (8.5) 

The above result can be rewritten as an exact solution of QED for the positronium ground 
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state in the nonrelativistic limit. If we let the color factor C/ = 1 and tuq rrie, g ^ e, 
eq. (|8.5|) can be reexpressed in terms of the Bohr radius and Bohr energy (in unit h = 1): 

«o = ^ = Eo = A = -— , « = ;r 8.6 

ZLJo niee^ 2ao An 

where the Bohr radius is redefined since we use a Gaussian-type trial wavefunction which is 
not the same as the exact hydrogen atomic ground state wavefunction. 

In fact, if we did not assume the A-independent of the QED coupling constant, we would 
obtain, except for a color factor, the same similarity renormalization group f3 function of 
eq. (|7.10|) . Thus, the fixed point of QED running coupling constant, e = 0, becomes infrared 



unstable, which is inconsistent with the well-known perturbative QED result of being an 
infrared stable fixed point. It implies that the similarity renormalization scheme can be 
applied to QED only for A = 0. With A = 0, our formulation reproduces the well-known 
method for solving QED bound states, namely, the nonperturbative bound states is deter- 
mined by solving the Schrodinger equation with the Coulomb potential (because the effective 
Hamiltonian Hxq only contains the Coulomb interaction at A = 0) and the remaining rela- 
tivistic radiative corrections (described here by Hxi) can then be systematically computed 
in perturbation theory. In other words, the limit A brings the effective theory back to 
the full theory in QED. 

However, unlike the QED, we cannot assume that the coupling constant g in QCD is 
independent of the scale A. In other words, one cannot freely take the limit of A — 0. 
With A being decreased, g\ becomes larger and larger. Thus, the resulting Hxo containing 
the Coulomb interaction alone is not sufficient to describe QCD bound states since the cor- 
rections from Hxj cannot be handled perturbatively. Therefore we are in practice unable 
to write down the full QCD theory in the weak-coupling formulation. The weak-coupling 
treatment of QCD in the limit A ^ is no longer valid. Indeed, we find that with a small 
A (A < 0.4 GeV), eq.( |6.23| ) has no solution. In other words, the so-called full QCD theory 
with massless gluon in low energy domain may only be formally interesting. To reproduce 
the acceptable hadronic properties with a small A value, we must include more complicated 
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higher order contributions from H\i into the nonperturbative bound state equation. Since 
the lowest bound value of A is the constituent gluon mass, a finite A, namely a nonzero con- 
stituent gluon mass, has effectively moved the nonperturbative contribution in the higher 
order processes into the low energy two-body confining interaction. The fact that the con- 
fining interaction dominates the binding dynamics of the quarkonium bound states at a 
finite A (about 1 GeV) is indeed an evidence why nonperturbative QCD can be treated as 
a weak-couple problem in the present formulation. The limit A — that can bring the 
effective theory back to the full theory is only implicative in QCD. 

On the other hand, the result from quarkonium ground states seems to tell us that with 
the larger A, the smaller the effective coupling constant can be reached. But this does not 
imply that the weak-coupling treatment to nonperturbative QCD works better for a larger 
A. The scale dependence of the wavef unction provides a restriction on the value of A. For 
quarkonia, uj\ is decreased with increasing A. However, uj\ is proportional to the mean 
value of the (transverse) momentum square of the quark inside the heavy quarkonia which 
characterizes the size of hadrons. Therefore it should not be too large or too small in the 
best description for bound states. For the range of eq.(|8.1|), we have. 



uxQ = 0.24 ~ 0.2 GeV. (8.7) 

Correspondingly, 

(r) ~ — = 0.8 ~ 1.0 fm, (8.8) 

which gives a resonable quarkonium size. Therefore, the true low energy QCD theory is 
determined by Hx with A being around the hadronic mass scale. 

Although our formulation thus explicitly involves A, by the requirement of the similarity 
renormalization group invariance, all the physical observables computed in this effective 
theory can still be A independent. The naive limit of A ^ that brings the effective theory 
back to the full QCD theory may imply that the final physical results calculated from H\ 
could also be gauge and Lorentz invariant, although Hx itself does not have these symmetries 
(consequently the hadronic wavefunctions must also be A-dependent). 

61 



In conclusion, the weak-coupling treatment approach to hadronic bound states in QCD 
can work well with the cutoff A being around the hadronic mass scale. The well-defined 
bound state approach in QED is a special case (A = 0) of the similarity renormalization 
group approach. The whole idea of the weak-coupling treatment to nonperturbative QCD 
via the similarity renormalization group approach is originally motivated from the bound 
state description of QED [|l|. Now, a consistent connection between QCD and QED and 
their differences in low energy domain is explicitly examined. 

IX. SUMMARY 

In this last section, we shall summarize the general formulation and the main results we 
have obtained, and then briefly discuss the further works. 

To realize the weak-coupling treatment of nonperturbative QCD recently proposed by 
Wilson and his collaborators |1|], we have studied explicitly the heavy quark bound state 
problem, based on the recently developed light-front heavy quark effective theory of QCD 
I^J^. Firstly, we have used the similarity renormalization group approach to derive the 
effective confining Hamiltonian in the low energy scale for heavy quarks in heavy mass limit. 
To make the similarity renormalization approach practically manable, we have introduced a 
local cutoff scheme (|2.4| ) to the bare QCD (and the effective heavy quark) Hamiltonian, which 
simplifies the cutoff scheme in [^. Meanwhile we have also introduced a simple smearing 
function fxij ( p.6|) to the similarity renormalization group approach which further simplifies 
the original formulation of The resulting low-energy effective QCD Hamiltonian of 
heavy quark interactions exhibits the coexistence of a confining interaction and a Coulomb 
interaction on the light-front. 

The realization of the weak-coupling treatment to nonperturbative QCD dynamics is 
very much based on the reseparation of the low energy effective Hamiltonian Hx into a 
nonperturbative part Hxo and the remaining as a perturbative term Hxi, and also on the 
use of the constituent picture, as we have seen throughout the present work. The use of the 
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constituent picture in light-front field theory allows us to expand the heavy hadrons only 
with the valence quark Fock space. The light-front heavy quark effective theory also largely 
simplifies the structure of the heavy hadron bound state equations [see ( [4.26| ) and ( |4.31| )]. 



A true realization of the weak-coupling approach to nonperturbative QCD can be ob- 
tained after solving the light-front bound state equations, from which one can determine 
the A-dependence of the effective coupling constant in H\, as a solution of the similarity 
renormalization group invariance. We have used a well-behaved light-front wavefunction 
ansatz (a Gaussian form) to analytically solve the quarkonium bound state equation and 
determine the scale-dependence of the effective coupling constant. We have also shown that 
the effective coupling constant at the hadronic scale A can be arbitrarily small. Thus, the 
possible weak-coupling treatment to nonperturbative QCD proposed in ref. is explicitly 
achieved. 

The results obtained in this paper is very optimistic. First, the A-dependence of the 
effective coupling constant determines the similarity renormalization group (3 function, from 
which some qualitative running behavior of the coupling constant in low energy domain 
is obtained. The running coupling constant ( |7.16D is qualitatively the same one used in 



the successful Richardson QQ potential for small momentum transfer (large distance). The 
later is a basic assumption to ensure the existence of a linear confining QQ potential in 
large distance, which is now obtained from QCD in the present work. A light-front picture 
of quark confinement from QCD is naturally manifested in iJ^o- If encompasses the general 
properties of these phenomenological confining potentials, the Richardson potential, the 
linear plus Coulomb potential and the logarithmic potentials used in the phenomenological 
description of quarkonium dynamics. 

The weak-coupling treatment can be realized for nonperturbative QCD because the simi- 
larity renormalization group approach with a finite A has extracted the confining interaction 
from the higher order nontrivial quark-gluon interactions into i^Ao- Equivalently speaking, 
the similarity renormalization group approach implicatively provides an inverse strong in- 
teraction expansion of QCD via the low energy cutoff scale A. As a physical consequence, 
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the confining interaction plays a dominant role in hadronic bound states, as we have seen in 
Section VII. The weak-coupling treatment of nonperturbative QCD is manable for A being 
around the hadronic mass scale. The similarity renormalization group invariance can remove 
the A-dependence in all the physical observables obtained from the effective H\. The possi- 
ble connection of the effective theory to the full QCD theory in the limit A ^ may further 
imply that the final physical results obtained from Hx could be gauge and Lorentz invariant. 
Meanwhile, we have also shown that in consistence with the behavior of the QED (3 function, 
the similarity renormalization approach can be applied to QED only if A = 0. Therefore, 
the confining picture obtained for QCD does not exist in QED but the well-known QED 
bound state method is reproduced. Thus, the consistency of the weak-coupling formulation 
has been qualitatively examined. 

The applications of the present theory to heavy hadron spectroscopy and various heavy 
hadron decay processes can be simply achieved by numerically solving the bound state 
equations ( [4.26| ) and ( ^.31| ), and by further including the l/wig corrections (which naturally 
leads to the spin splitting interactions). These will be presented in a forthcoming paper pS 



Finally, we should analyze the systematical approximations used in the whole computations 
in this paper, and then discuss the further works along this direction. 

The entire derivations presented in this paper are purely based on the first principle 
of QCD. Applying to heavy hadrons, we took the heavy mass limit so that the QCD is 
reduced to HQET but the first order kinetic energy has been also included in the leading 



order Hamiltonian for heavy quarkonia. In the forthcoming paper [^, the l/mq correc- 
tions (which contains all the spin splitting interactions) will be considered in the effective 
Hamiltonian Hxj in order to compute the heavy quarkonium spectroscopy and the next 
order correction to the bound states. These corrections should not affect on the main con- 
clusions obtained in this paper since the small factor 1 / itlq plus the weak coupling constant 
guarantee that these l/mg corrections can be computed perturbatively with respect to the 
nonperturbative heavy hadron bound states. 

The low energy QCD Hamiltonian H\ is obtained by a similarity renormalization trans- 
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formation to the bare QCD Hamiltonian (bare HQET Hamiltonian for heavy quarks). With 
the idea of the weak-couphng treatment to low energy QCD dynamics, the nonperturbative 
part if AO is computed upon to the second order in H\ and the bound states are truncated 
to only including the valence quark Fock space. The higher order corrections in Hx (in- 
cluded both the l/rriQ corrections and the radiative corrections) can be examined in the 
usual Hamiltonian perturbation theory which will also provide a consistent check to 
the validity of the present weak-coupling treatment. With these corrections computed in 
the Hamiltonian formulation, the contributions from the higher Fock space will be naturally 
included. 

Finally, instead of numerically solving the light-front bound state equations, we used here 
a well-behaved light-front wavefunction ansatz to determine the bound state equation for 
heavy quarkonia. A numerical calculation of the light-front bound state equations ([4.26|) and 
( [4.31| ) for the heavy hadronic spectroscopy will be also presented in ||2^. Such a numerical 
computation of the bound states is much simpler in comparison to the lattice QCD simulation 

and it will also directly give the hadronic wavefunctions in physical space. The resulting 
wavefunctions of the quarkonium bound states combining with the systematic computations 
of the subsequent radiative and l/mg corrections thus provide a truly unified first- principles 
QCD description for various heavy quarkonium annihilation and production processes. 

All the derivations presented in this paper are rigorous QCD derivation in low energy 
domain. The extension of the computations to heavy-light quark systems is straightforward 
but is more attractive in the current investigations on heavy hadrons. The extension of the 
present work to light-light hadrons requires a understanding of the chiral symmetry breaking 
mechanism in QCD which is a new challenge to nonperturbative QCD on the light-front. 
Nevertheless, the present work has provided a preliminary realization to the weak-coupling 
treatment of nonperturbative QCD proposed recently by Wilson et al. [|I| . The new research 
along this direction is now in progress. 
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FIGURES 

FIG. 1. A plot of the confining potential energy, the Coulomb potential energy and the total 
potential energy as functions of the dimensionless variable xu but w is proportional to ~ ri via ijJ\. 

„2y2 

The energies are scaled by the factor -^;^Cf. 

FIG. 2. The A-dependence of the effective coupling constant a\. The solid line is given by 
the analytical result (^^), and the dots are obtained by numerically minimizing the quarkonium 
binding energy ( 6.231 ). Here A is given in units of GeV. 



FIG. 3. The A-dependence of the wavefunction parameter ujxq which is the solution of mini- 
mizing the quarkonium binding energy (|6.23|), where A is given in units of GeV. 
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TABLES 



Table I. Solution for charmonium ground state with rric — lA GeV 



A (GeV) 


ax 


00X0 (GeV) 




0.2 


0.02665 


0.222 


0.026836 


0.053173 


0.080009 


0.3 


0.06480 


0.275 


0.067902 


0.112018 


0.179920 


0.4 


0.11831 


0.314 


0.130225 


0.189781 


0.320006 


Tal)le II. Solution for l)ottomonium ground state with riii, = 1.8 GeV 


A (GeV) 


ax 


^Ao (GeV) 




0.15 


0.029965 


0.492 


0.023397 


0.021602 


0.044999 


0.20 


0.06795 


0.623 


0.050183 


0.029816 


0.079999 


0.25 


0.1385 


0.779 


0.098074 


0.026930 


0.125004 



Table III. Some numerical solution on the A-dependence 



of the weak coupling constant ax- 





charmonium 


bottomonium 


A (GeV) 


ax 


00x0 (GeV) 


ax 


00x0 (GeV) 


0.75 


0.05960 


0.241 






1.0 


0.02665 


0.222 


0.06795 


0.623 


1.5 


0.00912 


0.199 


0.01607 


0.478 


2.0 


0.00441 


0.185 


0.00695 


0.427 
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Table IV. The A-dependence of various interactions to the binding energy 



A (GeV) 


0.5 


0.75 


1.0 


1.2 


1.4 


1.8 


2.0 


3.0 


Skin (GeV^) 


0.04 


0.031 


0.027 


0.025 


0.023 


0.021 


0.02 


0.018 


Sconf (GeV^) 


0.049 


0.050 


0.053 


0.055 


0.057 


0.059 


0.06 


0.062 


Scoi (GeV2) 


-0.009 


-0.001 


-0.001 


0.0 


0.0 


0.0 


0.0 


0.0 
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Fig.2 





Fig.l 



